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ABSTRACT 


Current  two-dimensional  preliminary  design  codes  use  structured  programming, 
which  is  rigid  and  does  not  allow  the  user  to  vary  parameters  easily.  This  study  uses 
object-oriented  programming  to  allow  the  user  to  vary  all  selectable  parameters  in  a 
familiar  Windows  operating  environment.  The  programmed  design  is  based  on  the 
assumptions  of  axial  and  free- vortex  flow  between  blade  rows,  simple  radial  equilibrium, 
and  a  thermally  and  calorically  perfect  gas.  The  program  allows  a  fan  or  core  stage 
design  and  uses  an  open  architecture  to  facilitate  upgrades  and  extensions. 

Using  the  Naval  Postgraduate  School’s  (NPS)  transonic  compressor  design  as 
input,  the  preliminary  design  code  output  was  compared  to  the  detailed  throughflow 
design  of  the  transonic  compressor.  The  results  agreed  reasonably  well  with  detailed 
throughflow  design.  With  some  minor  improvements  this  code  can  easily  be  used  to 
develop  a  preliminary  design  that  can  be  optimized  to  the  user’s  requirements. 
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1.  INTRODUCTION 


The  design  of  a  new  axial  compressor  involves  a  sequence  of  steps,  progressing 
through  a  sequence  of  computational  programs  of  increasing  complexity  and 
sophistication.  The  first,  or  ‘preliminary  design’  step,  can  be  a  one -dimensional  ‘mean¬ 
line  calculation’,  or  a  two-dimensional  calculation  of  a  preliminary  flow  path  and 
selection  of  the  blading.  The  latter  is  the  minimum  required  if  the  overall  task  is  the 
preliminary  design  of  an  aircraft  gas  turbine  engine.  It  is  also  what  is  required  to  provide 
the  input  for  detailed  throughflow  design  codes  [Ref.  1],  which,  in  turn,  generate  inputs  to 
codes,  which  compose  the  blade  geometry  [Ref.  2]  for  manufacturing. 

The  current  preliminary  compressor  design  code  used  in  aircraft  engine  design 
courses  [Ref.  3],  was  developed  progressively  using  different  versions  of  Hewlett- 
Packard  (HP)  BASIC  [Ref  4] .  This  highly  structured  programming  language  is  rigid  and 
does  not  allow  the  user  to  vary  one,  or  several,  parameters  easily  in  order  to  change  or 
optimize  a  design.  Use  of  the  programs  requires  the  installation  of  HP  BASIC  for 
Windows. 

The  purpose  of  the  present  study  was  to  develop  a  preliminary  compressor  design 
code  that  would  satisfy  the  following  conditions: 

•  Be  simple  to  use. 

•  Allow  all  selectable  parameters  to  be  changed. 

•  Require  only  the  Windows  operating  system. 

•  Use  an  open  architecture  to  allow  upgrades  (e.g.,  different  loss  models)  or 
additions  (e.g.,  turbine  design) 

To  meet  these  conditions,  Microsoft’s  Visual  Basic  6.0  [Ref.  5]  was  selected  and 
used  to  develop  the  source  code. 
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II.  PROGRAM  OVERVIEW 


A.  ASSUMPTIONS 

The  following  assumptions  were  made  in  the  development  of  the  programmed 
equations  and  design  of  the  code: 

•  Design  is  for  an  axial  flow  compressor  (with  radial  movement  of  the  mean 
line). 

•  Simple  radial  equilibrium  is  assumed  from  hub  to  tip. 

•  Free  vortex  flow  is  assumed  between  blade  rows. 

•  A  conceptual  engine  design  study  [using  Ref.  6  for  example]  will  generate  the 
inputs  for  the  code. 


B.  INPUTS  AND  OUTPUTS 

The  stage  is  shown  schematically  in  Fig.  1,  and  the  inputs  required  by  the  code, 
are  shown  in  Table  1. 
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Inlet  (1) 

Rotor  (2) 

Stator  (3) 

m 

D2m 

A31 

Pt 

a  2m 

a  3m 

Tt 

R2I 

R32 

R 

<1^21 

<1^32 

Y 

52 

63 

0) 

AR2 

AR3 

Mzit 

t/C2h 

t/C3h 

Mwit 

t/C2m 

t/C3ni 

ait 

t/C2t 

t/C3t 

fo2 

fa3 

Z2 

Z3 

Tal 

ble  1.  User- Provided  Inputs 

At  the  inlet  station,  it  is  required  to  specify  the  gas  (through  the  gas  constant  and 
ratio  of  specific  heats;  a  thermally  and  calorically  perfect  gas  is  assumed),  the  mass  flow 
rate,  and  the  stagnation  conditions.  Then,  in  order  to  accommodate  the  usual  fan  or  core 
design  constraints,  four  parameters  must  be  specified  for  the  rotor;  namely,  rotational 
speed,  axial  Mach  number,  relative  Mach  number,  and  flow  angle  at  the  tip.  Only  three 
of  these  four  are  independent.  Off-  line  calculation  is  required  to  determine  the  other  one. 

The  design  selections  for  the  rotor  and  stator  are  listed  in  columns  two  and  three 
of  Table  1.  Note  that  the  blade  loading  throughout  the  stage  is  determined  by  the 
selection  of  the  diffusion  factor  at  only  one  location  for  one  blade.  Blade  aspect  ratio  and 
thickness  variations  are  choices,  which  are  determined  by  structural  considerations. 
Structural  constraints  are  not  included  in  the  program. 

The  parameters,  which  are  calculated  and  output  by  the  code,  are  shown  in 
Table  2. 
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Inlet  (1) 

Rotor  (2] 

Stator  (3] 

Comp. 

Hub 

Mean 

Tip 

Hub 

Mean 

Tip 

Hub 

Mean 

Tip 

(h) 

(m) 

(t) 

(h) 

(m) 

(t) 

(h) 

(m) 

(t) 

X 

Xih 

Xim 

Xit 

X2h 

X2m 

X2t 

X3h 

X3m 

X3t 

Xz 

Xzih 

Xzim 

Xzit 

Xz2h 

Xz2m 

Xz2t 

Xz3h 

Xz3m 

Xz3t 

Xu 

Xuih 

Xulm 

Xuit 

Xu2h 

Xu2m 

Xu2t 

Xu3h 

Xu3m 

Xu3t 

Xw 

Xwih 

Xwim 

Xwit 

Xw2h 

Xw2m 

Xw2t 

Xw3h 

Xw3m 

Xw3t 

Xe 

Xeih 

Xeim 

Xeit 

Xe2h 

Xe2ix 

Xe2t 

Xe3h 

Xe3m 

Xe3t 

r 

rih 

rim 

rit 

r2h 

r2m 

r2t 

r3h 

r3m 

r3t 

M 

Mih 

Mim 

Ml, 

M2h 

M2m 

M2t 

M3h 

M3m 

M3t 

Mz 

Mzih 

Mzim 

Mzit 

Mz2h 

Mz2m 

Mz2t 

Mz3h 

Mz3m 

Mz3t 

Mw 

Mwih 

Mwim 

Mwit 

Mw2h 

Mw2m 

Mw2t 

Mw3h 

Mw3m 

Mw3t 

P 

Plh 

Plm 

Pit 

P2h 

P2m 

P2t 

P3h 

P3m 

P3t 

a 

ttih 

OClm 

Input 

0(2h 

0(2  m 

a2t 

0(3h 

0(3  m 

a3t 

Y 

Y2h 

Y2m 

Y2t 

Y3h 

Y3m 

Y3t 

Yw 

Yw2h 

Yw2m 

Yw2t 

Yw3h 

Yw3m 

Yw3t 

D 

D2h 

Input 

D2t 

D3h 

D3m 

D3t 

o 

02h 

Input 

02t 

(53h 

Input 

(53t 

Tst 

rst2h 

rst2m 

rst2t 

rst3h 

rst3m 

rst3t 

1* 

t*2h 

l*2m 

l*2t 

l*3h 

l*3m 

l*3t 

(|)*2h 

(^*2m 

4**21 

4)*3h 

4>*3m 

4>*3t 

5* 

5*2h 

5*2m 

5*2t 

5*3h 

5*3m 

5*3t 

rht 

rhti 

rht2 

rht3 

A 

Ai 

A2 

A3 

T/Tti 

Ti/Tti 

T2/r,i 

T3/rti 

P/Pti 

Pi/Pti 

P2/Pt2 

P3/PG 

Tt/Tti 

Tt2/Tti 

Tt3/rti 

(X) 

Pt/Pti 

Pt2/Ptl 

Pt3/Ptl 

(n) 

A/Ai 

A2/A1 

A3/A1 

0)p 

(0p2 

(Op  3 

(Osftc 

(0sftc2 

(0sftc3 

(Os 

C0s2 

(0s3 

cot 

(0x2 

(0x3 

$ 

^2m 

^3m 

H 

H2 

H3 

C 

C2 

C3 

S 

S2 

S3 

^ev 

Zrev2 

Zrev3 

ARrev 

ARrev2 

ARrev3 

Crev 

Crev2 

Crev3 

Table  2.  Outputs 
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m.  ALGORITHMS 


A.  PROGRAM  SEQUENCE 

The  overall  program  sequence  is  shown  in  Fig.  2.  After  the  program  is  started  a 
“splash”  screen  is  first  shown.  A  splash  screen  is  an  introductory  screen  (similar  to  an 
“about”  screen)  which  states  the  name,  owner(s),  and  version  of  the  program.  Next,  an 
interactive  screen  appears,  giving  the  user  a  choice  of  either  going  through  a  compressor 
design  or  a  turbine  design.  The  scope  of  this  paper  covers  only  a  compressor  design; 
therefore  the  turbine  option  has  been  disabled.  When  the  Compressor  Design  button  is 
pressed  the  input  screen  will  appear.  The  design  sequence  is  shown  in  Fig.  3.  The 
required  inputs  are  shown  in  Table  1.  Once  the  inputs  are  typed  in  the  user  presses  the 
OK  button  and  the  main  screen  appears.  The  main  screen  is  laid  out  in  a  tab  format  with 
five  tabs.  The  tabs  are  as  follows: 

•  Inlet  Conditions  (1) 

•  Rotor  Calculations  (2) 

•  Stator  Calculations  (3) 

•  Stage  Performance 

•  Blade  Geometry 

This  allows  the  variables  and  code  to  be  grouped  so  the  user  can  better  understand 
what  is  being  displayed.  It  also  allows  the  programmer  to  develop  the  code  with  an  open 
architecture  for  easy  updates.  At  this  point  the  user  can  still  modify  the  inputs  before 
executing  the  design  calculations. 
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Figure  2.  Main  Program  Flowchart 
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Design  Flowehart 
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The  user  has  the  option  to  either  apply  the  inputs  to  the  design  calculations  or 
calculate  the  next  stage.  A  multistage  approach  is  also  outside  the  scope  of  this  paper  so 
this  option  has  been  disabled.  After  the  user  presses  the  Apply  button  the  program 
calculates  the  basic  throughflow  velocity  diagrams,  losses,  performance  and  geometry. 
During  the  throughflow  calculations  an  input  screen  is  displayed  for  the  user  to  input  the 
whole  number  of  blades  needed  (Z)  for  the  rotor  and  stator.  The  main  screen  is  displayed 
again  with  the  outputs  filled  in  the  text  boxes.  The  outputs  are  shown  in  Table  2.  With 
the  initial  throughflow  calculations  complete  the  user  may  end  the  program  of  modify  the 
inputs  until  the  desired  results  are  achieved. 

B.  DESIGN  EQUATIONS 

The  equations  programmed  in  the  source  code  are  listed  in  Appendix  A.  The 
equations  are  grouped  the  same  as  the  tabs  on  the  main  screen.  The  equations  for  the 
velocity  diagrams,  profile  loss,  secondary  flow  loss,  tip  clearance  loss,  and  stage 
performance  are  from  Shreeve  [Ref.  7  and  8].  The  shock  loss  equation  is  from  Koch  and 
Smith  [Ref.  9]. 

The  incidence,  deviation  and  camber  angles  (which  relate  the  flow  angles  to  the 
blade  geometry,  as  shown  in  Fig.  4)  were  derived  from  NASA  SP-36  [Ref.  10].  Sixth 
degree  polynomial  curvefits  were  used  to  approximate  the  data  in  Figs.  137,  138,  142, 
161,  162,  168,  172,  178,  179,  and  180  of  Ref.  10.  An  interpolation  routine  [Ref.  11],  was 
programmed  to  solve  for  unknowns  within  the  curve- fits. 
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The  performance  of  the  stage  is  calculated  following  the  thermodynamic  process 
shown  in  Fig.  5.  The  conditions  on  the  mean  line  are  taken  as  being  representative  of  the 
stage;  however,  the  loss  coefficients  include  contributions  due  to  secondary  flow  and  tip- 
clearance. 
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Figure  5.  T-s  Diagram  for  a  Compressor  Stage 
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IV.  PROGRAM  STRUCTURE 


A.  OBJECTS 

Forms  (a.k.a.  screens)  are  used  to  interact  with  the  user.  Forms  are  ‘objects’  in 
Visual  Basic  and  Visual  C++  programming.  Other  common  objects  include  text  boxes, 
option  boxes  and  tabs.  Objects  have  properties  (this  is  how  the  object  looks  to  the  user) 
which  can  be  affected  during  design  or  run  time  (program  execution).  The  objects  allow 
the  user  to  modify  the  inputs  throughout  the  program  and  either  executes  the  desired 
modification  or  keeps  the  original  inputs.  The  objects  (or  “screen  captures”)  are  shown 
in  Appendix  B.l. 

B.  VARIABLES 

Variables  are  used  in  programming  as  placeholders  of  data  where  they  are  used 
throughout  the  code  to  execute  statements.  Two-dimensional  arrays  are  used  extensively 
throughout  the  code  since  they  allow  a  parameter  that  changes  both  throughflow  and 
spanwise  to  be  stored  in  one  variable.  For  example,  (relative  velocity)  varies  from 
inlet  (1)  to  stator  (3)  and  from  hub  (1)  to  tip  (3).  If  we  store  this  variable  separately  we 
would  need  to  track  nine  variables.  Instead,  we  simply  track  one  variable  Xw(3,3)  where 
the  first  number  is  the  throughflow  number  and  the  second  number  is  the  spanwise 
number.  Appendix  B.2  shows  all  the  variables  used  in  the  program.  This  includes  all  the 
input  and  output  parameters  as  well  as  others  that  are  not  seen  by  the  user. 


C.  MODULES 

Code  is  used  to  state  what  needs  to  be  executed.  The  code  syntax  is  based  on 
BASIC,  which  translates  “pseudo”  language  to  machine  language.  Code  is  written  for 
objects  to  tell  the  object  what  to  do  when  acted  upon  by  the  user.  For  this  program  there 
is  very  little  code  for  objects.  Most  of  the  code  is  a  sequence  of  design  equations  which 
do  not  cause  changes  in  any  objects  and  require  minimal  interaction  with  the  user  during 
execution.  Modules  are  used  in  this  situation.  Modules  take  subroutines  or  functions  that 
do  not  affect  an  object  and  keep  them  in  a  separate  file  for  easier  reading  during 
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programming  or  debugging.  The  code,  under  the  form  ‘frmCompressor’,  is  the  design 
code  sequence  shown  in  Figure  2.  Appendix  B.3  shows  all  the  source  code  for  the 
program. 
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V.  RESULTS  AND  DISCUSSION 


Sanger’s  transonic  compressor  design  [Ref.  13]  was  used  to  test  the  ability  of  the 
code  to  approximate  a  known  axial  stage  design.  A  set  of  hand  calculations  was  also 
carried  out  independently  of  the  code  in  order  to  both  validate  the  coding,  and  to 
document  the  test  case.  The  design  input  flow  conditions,  and  parameters  derived  from 
the  final  geometry  of  the  Sanger  design,  were  used  as  inputs  to  both  the  preliminary 
design  code  and  hand  calculations.  Appendix  C.l  and  Appendix  C.2  document,  in  detail, 
the  results  of  the  hand  calculations  and  the  preliminary  design  code,  respectively. 

The  results  from  hand  calculations  and  from  the  code  were  compared  to  the 
output  of  the  streamline  curvature  code  applied  to  the  Sanger  design,  which  is  given  in 
Appendix  E  of  Sanger  [Ref.  13].  The  comparisons  are  shown  in  Appendix  D.  It  can  be 
seen  from  the  comparison  charts  in  Appendix  D  that  the  hand  calculations  agreed  fully 
with  the  preliminary  design  code  calculations.  Also,  for  most  parameters,  the  hand  and 
preliminary  design  code  calculations  agreed  with  the  streamline  curvature  code  outputs. 
In  Figure  Dl,  the  calculated  annulus  geometry  agrees  well  at  the  inlet  and  then  begins  to 
deviate  somewhat  through  the  stage.  This  is  because  the  streamline  curvature  code  takes 
into  account  blockage,  whereas  the  preliminary  design  code  and  hand  calculations,  do 
not.  This  is  easy  to  correct.  The  stage  performance  however,  shown  in  Fig.  D2,  is 
predicted  very  successfully  by  the  code.  In  Fig.  D3  and  D4,  the  differences  in  solidity  and 
blade  height  are  also,  indirectly,  the  result  of  omitting  blockage  from  the  calculation  of 
annulus  area.  Velocity  diagram  details  are  compared  in  Figs.  D5  to  D8.  It  is  clear  that  the 
preliminary  design  code  reproduces  the  final  design  values  to  very  acceptable  accuracy. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


From  the  comparisons  given  in  Appendix  D,  it  can  be  seen  that  the  preliminary 
design  code  does  very  well  in  developing  the  velocity  diagrams  and  the  initial  blading 
geometry  necessary  for  a  detailed  throughflow  design  and  final  geometry  calculation 
[Ref.  3  and  4]. 

Improvements  can  be  made  in  order  to  have  the  design  more  detailed  as  well  as 
improve  the  code’s  usability.  They  are  as  follows: 

•  Add  the  ability  to  do  a  multi-stage  design  (use  the  output  of  the  previous  stage 
as  an  input  to  the  new  stage). 

•  Add  different  input  screens  for  specific  design  cases  (e.g.,  fan  or  core). 

•  Add  stress  limits  (hoop  and  centrifugal)  for  fan  design. 

•  Add  blockage  and  bleed  to  the  throughflow  calculations. 

•  Draw  scale  velocity  diagrams. 

•  Draw  blades  based  on  code  output  (i.e.,  built  in  geometry  package). 

•  Add  parametric  analysis  for  a  range  of  values  (e.g.,  inlet  flow  angle). 

•  Add  turbine  stage  design  (in  parallel). 

•  Add  the  ability  to  open  and  save  data. 

•  Add  the  ability  to  print  the  user’s  results. 

•  Compile  the  code  into  a  stand-alone  executable. 

Incorporation  of  these  improvements  will  make  the  code  a  preliminary 
turbomachine  design  software  package  that  can  be  used  as  inputs  to  detailed  design 
packages,  as  well  as  providing  a  needed  teaching  tool  for  aircraft  engine  design. 
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APPENDIX  A.  DESIGN  EQUATIONS 

A.l  shows  the  throughflow  equations  and  A. 2  shows  the  interpolation  equations 
for  the  incidence,  camber  and  deviation  angles. 


THROUGHFLOW  EQUATIONS 
1.  Inlet  Conditions 

Given:  m  ,  Pt,  Tt,  Mwit>  CO ,  R,  y 
Vary:  a.,  (0<a„  <P.,) 
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Due  to  Radial  Equilibriums 
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Initially  set  ,  cbj  =  0 
On  subsequent  iterations  cbsprc,  and  are 
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5, 


Stator  Performance 
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A3  _  4>[,  cosa,,  Vx 
A,  03„-cosa3„  j 
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I-H3 
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6.  Rotor  &  Stator  at  the  Hub  and  Tip 

The  below  equations  are  generic  forms  of  the  specific  equation.  The  brackets  []  replace  the  rotor 
(2)  and  stator  (3)  subscripts.  Also,  the  parentheses  ()  replace  the  hub  (1)  and  tip  (3)  subscripts. 
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Mno=' 


Y-1 


1-y: 


Mwno  =' 
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-YJ, 


1-Y 

K  wno 

Blade  Geometry 


The  below  equations  are  generic  forms  of  the  specific  equation.  The  brackets  []  replace  the  rotor 
(2)  and  stator  (3)  subscripts. 
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8.  Efficiency 
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A.2.  INTERPOLATION  EQUATIONS 

Starting  with  a  general  quadratic  equation 
D‘(R)  =  Ai  R"+  R+q 


solve  for  Aj,  Bj,  and  Q  using  known  points  i-1,  i,  and  i+1. 
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Over  the  first  and  last  intervals,  only  one  quadratic  can  be  defined,  so  that 

D, „(i)=Y-(R:-R:)+Y-(R:-R:)+q  ■(^ 

D,„.  (N)  =  (RL  -R;)+  ^  (K, -K)+C,  (R,,.  - R, ) 


The  complete  integral  is  given  by 

Rn+1  N-1 

I  D(R)dr  =  D.„,(l)  +  J^D.„.(i)  +  D,.(N) 

R  i=2 


27 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


28 


APPENDIX  B.  PROGRAM  SOURCE  DATA 

B.l.  SCREEN  SNAPSHOTS 


Figure  Bl.  Splash  Screen 
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Figure  B3.  Initial  Inputs  Screen 


8a>  ^P*y  tW**  Ml 


- 3 


Cto  M 

Afipty  N$xt9tMf  ' 


m 


••HCiMWMuiWtti 


il’> 


9flC*  l7i 


Knowns 

m 

|1/(B 

p, 

|MS9 

A>:_ 

jiira 

Ceostanls  and 

f  Param*t«<« 

R: 

v; 

b4e?54 

€t.,  ' 

[i5 

|:S3B13 


Inlet  VetocHy  Diagram  - 


Re&UK 


I - 

a;. 

I 

A'.,.:l 

'A„:J 

A.:l 

I - 

A'™; 

1 - 

- 

Lv^ :  1 

A  1 

a;. 

r 

a;„  :  1 

'A'„  .  1 

«C  1 

I - 

A- 

r 

A„:l 

A/.;  1 

"»-i 

I 

■[ 

A'..:i 

A/,„ ;  i 

A',.1 

I 

: 

1 - 

A  „ :  i 

A/..:  J 

A',„i 

I 

r,  ■ 

1 

AC.:I 

AC.  i 

-C  I 

I - 

r,  '■ 

I 

A/.-.  1 

AC„I 

Figure  B4.  Compressor  Design  Screen:  Inlet  Conditions  (1) 
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Figure  B5.  Compressor  Design  Screen:  Rotor  Calculations  (2) 
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Figure  B6.  Compressor  Design  Screen:  Stator  Calculations  (3) 
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Figure  B7.  Compressor  Design  Screen:  Stage  Performance 


Figure  B8.  Compressor  Design  Screen:  Blade  Geometry 
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Figure  B9.  Blade  Number  Input  Screen 
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B.2.  KEY  VARIABLES 


Error!  Not  a  valid  link. 

Table  Bl.  Velocity  Diagram  Variables 


Error!  Not  a  valid  link. 

Table  B2.  Stage  Performance  Variables 


Error!  Not  a  valid  link. 

Table  B3.  Blade  Geometry  Variables 

B.3.  SOURCE  CODE 

Lines  with  an  apostrophe  (')  in  the  front  is  a  comment  line  in  the  code  and  is  not 
an  executed  statement. 

1.  Splash  Screen  Code 

Option  Explicit 
Private  Sub  Form_Load() 

IblVersion. Caption  =  "Version  ”  &  App. Major  &  &  App. Minor  &  &  App. Revision 

IblProductName. Caption  =  App. Title 
IblCompany. Caption  =  App.CompanyName 

End  Sub 

Private  Sub  tmrSplash_Timer() 

Unload  Me 
frmCompT  urb  .Show 

End  Sub 

2.  Compressor/Turbine  Selection  Screen  Code 

Private  Sub  cmdComp_Click() 

Unload  Me 
frmlnitial.Show 

End  Sub 

3.  Initial  Input  Screen  Code 

Option  Explicit 
Option  Base  1 

Private  Sub  OKButton_Click() 
frmlnitial.Hide 
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frmCompressor.Show 
'  Transfer  initial  values  to  the  main  form. 

End  Sub 

4.  Compressor  Design  Screen  Code 

Option  Explicit 
Option  Base  1 
Public  i  As  Integer 

'  Integer  to  use  as  an  array  counter  from  inlet  (i  =  1)  to  rotor  (i=2)  to  stator  (i=3) 

'  It  should  be  noted  that  although  some  arrays  have  no  inlet  values  (e.g.  diffusion  and  deg.  of 
reation) 

'  we  still  set  the  array  value  to  3  versus  2.  This  is  to  keep  consistent  with  the  numbering. 

'  arrays  that  don't  have  the  i=l  value  assigned  are  null  and  irrelevent  since  they  are  not  used 
'  in  any  calculations. 

Public  j  As  Integer 

'  integer  to  use  as  an  array  counter  from  hub  (j  =  I)  to  mean  (j  =  2)  to  tip  (j  =  3) 

'  It  should  be  noted  that  although  hub,  mn  and  tip  have  been  set  up  a  constants 
'  j  will  be  used  to  calculate  the  minimum  loss  incidence  angle  and  the  camber  angle 

Public  ksh  As  Single 
'  correction  factor  for  shape 
Public  slopegraph  As  Boolean 

'  checks  which  graph  should  be  used  to  calculate  the  slope  factor 

Const  testvar  As  Single  =  0.0001 

'  exit  criteria  for  main  loop 

Public  loopcount  As  Integer 

'  loop  counter  to  track  the  number  of  iterations 

Dim  tempvar(3)  As  Double 

'  holder  for  previous  loss  value  for  comparison 

Public  CompEff  As  Double 

'  compressor  efficiency 

Private  Sub  cmdApply_Click() 

Sequence 
Update  VelDiag 
UpdatePerf 
UpdateGeo 

End  Sub 
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Private  Sub  Form_Initialize() 

frmCompressor.txtmdot.Text  =  frmlnitial.txtmdot.Text 
frmCompressor.txtPtl.Text  =  frmlnitial.txtPtl.Text 
frmCompressor.txtTtl.Text  =  frmlnitial.txtTtl.Text 
frmCompressor.txtR.Text  =  frmlnitial.txtR.Text 
frmCompressor.txtMzt(l).Text  =  frmlnitial.txtMzlt.Text 
frmCompressor.txtMwt(l).Text  =  frmlnitial.txtMwlt.Text 
frmCompressor.txtspeed.Text  =  frmlnitial.txtomega.Text 
frmCompressor.txtgamma.Text  =  frmlnitial.txtgamma.Text 
frmCompressor.txtalphat(l).Text  =  frmlnitial.txtAlphalt.Text 
frmCompressor.txtDm(2).Text  =  frmlnitial.txtDRm.Text 
frmCompressor.txtsigmam(2).Text  =  frmlnitial.txtsigmaRm.Text 
frmCompressor.txtsigmam(3).Text  =  frmlnitial.txtsigmaSm.Text 
frmCompressor.txtR21.Text  =  frmInitial.txtR21.Text 
frmCompressor.txtR32.Text  =  frmInitial.txtR32.Text 
frmCompressor.txtphi21.Text  =  frmInitial.txtphi21.Text 
frmCompressor.txtphi32.Text  =  frmInitial.txtphi32.Text 
frmCompressor.txtA31.Text  =  frmInitial.txtA31.Text 
frmCompressor.txtDelta(2).Text  =  frmlnitial.txtDeltar.Text 
frmCompressor.txtDelta(3).Text  =  frmlnitial.txtDeltas.Text 
frmCompressor.txtAR(2).Text  =  frmlnitial.txtARr.Text 
frmCompressor.txtAR(3).Text  =  frmlnitial.txtARs.Text 
frmCompressor.txttch(2).Text  =  frmlnitial.txttcrh.Text 
frmCompressor.txttch(3).Text  =  frmlnitial.txttcsh.Text 
frmCompressor.txttcm(2).Text  =  frmlnitial.txttcrm.Text 
frmCompressor.txttcm(3).Text  =  frmlnitial.txttcsm.Text 
frmCompressor.txttct(2).Text  =  frmlnitial.txttcrt.Text 
frmCompressor.txttct(3).Text  =  frmlnitial.txttcst.Text 
frmCompressor.txtfsigma(2).Text  =  frmlnitial.txtfsigmar.Text 
frmCompressor.txtfsigma(3).Text  =  frmlnitial.txtfsigmas.Text 
'  Initialize  Inlet  Condition  Variables. 
alpha(l,  3)  =  frmCompressor.txtalphat(l).Text 
gamma  =  frmCompressor.txtgamma.Text 
Mz(l,  3)  =  frmCompressor.txtMzt(l).Text 
mdot  =  frmCompressor.txtmdot.Text 
Mw(l,  3)  =  frmCompressor.txtMwt(l).Text 
omega  =  frmCompressor.txtspeed.Text 
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Pt(l)  =  frmCompressor.txtPtl.Text 
Rbar  =  frmCompressor.txtR.Text 
Tt(l)  =  frmCompressor.txtTtl.Text 
'  Initialize  Rotor  Assumption  Variables 
D(2,  2)  =  frmCompressor.txtDm(2).Text 
sigma(2,  2)  =  frmCompressor.txtsigmam(2).Text 
pitch(2)  =  frmCompressor.txtR21.Text 
phi21  =  frmCompressor.txtphi21.Text 
fsigma(2)  =  frmCompressor.txtfsigma(2).Text 
'  Initialize  Stator  Assumption  Variables 
A31  =  frmCompressor.txtA31.Text 
sigma(3,  2)  =  frmCompressor.txtsigmam(3).Text 
pitch(3)  =  frmCompressor.txtR32.Text 
phi32  =  frmCompressor.txtphi32.Text 
fsigma(3)  =  frmCompressor.txtfsigma(3).Text 
'  Initialize  Stage  Performance  variables 
delta(2)  =  frmCompressor.txtDelta(2).Text 
delta(3)  =  frmCompressor.txtDelta(3).Text 
'  initialize  blade  geometry  variables 
For  i  =  2  To  3 

For  j  =  1  To  3 

AR(i)  =  frmCompressor.txtAR(i).Text 
If  j  =  1  Then 

tc(i,  j)  =  frmCompressor.txttch(i).Text 
El  self  j  =  2  Then 

tc(i,  j)  =  frmCompressor.txttcm(i).Text 
El  self  j  =  3  Then 

tc(i,  j)  =  frmCompressor.txttct(i).Text 

End  If 

Next  j 

Next  i 

End  Sub 

Private  Sub  Update VelDiagO 
Eor  i  =  1  To  3 

Eor  j  =  1  To  3 

If  j  =  1  Then 

txtBetah(i).Text  =  beta(i,  j) 
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txtAlphah(i).Text  =  alpha(i,  j) 
txtXh(i).Text  =  x(i,  j) 
txtXUh(i).Text  =  Xu(i,  j) 
txtXzh(i).Text  =  Xz(i,  j) 
txtXthetah(i).Text  =  Xtheta(i,  j) 
txtXwh(i).Text  =  Xw(i,  j) 
txtrh(i).Text  =  R(i,  j) 
txtMh(i).Text  =  M(i,  j) 
txtMwh(i).Text  =  Mw(i,  j) 
txtMzh(i).Text  =  Mz(i,  j) 

If  i  =  2  Then 

txtrsth(i).Text  =  rst(i,  j) 
txtDh(i).Text  =  D(i,  j) 
txtYh(i).Text  =  Y(i,  j) 
txtYwh(i).Text  =  Yw(i,  j) 
El  self  i  =  3  Then 

txtrsth(i).Text  =  rst(i,  j) 
txtDh(i).Text  =  D(i,  j) 
txtYh(i).Text  =  Y(i,  j) 
txtYwh(i).Text  =  Yw(i,  j) 

End  If 

End  If 
If  j  =  2  Then 

txtbetam(i).Text  =  beta(i,  j) 
txtalphani(i).Text  =  alpha(i,  j) 
txtXm(i).Text  =  x(i,  j) 
txtXUm(i).Text  =  Xu(i,  j) 
txtXzm(i).Text  =  Xz(i,  j) 
txtXthetam(i).Tejt  =  Xtheta(i,  j) 
txLXwm(i).Text  =  Xw(i,  j) 
txtmi(i).Text  =  R(i,  j) 
txtMm(i).Text  =  M(i,  j) 
txtMwni(i).Text  =  Mw(i,  j) 
txtMzm(i).Text  =  Mz(i,  j) 

If  i  =  2  Then 

txtrstm(i).Text  =  rst(i,  j) 
txtYm(i).Text  =  Y(i,  j) 
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txtYwm(i).Text  =  Yw(i,  j) 

El  self  i  =  3  Then 

txtrstm(i).Text  =  rst(i,  j) 
txtYm(i).Text  =  Y(i,  j) 
txtYwm(i).Text  =  Yw(i,  j) 
txtDm(i)  =  D(i,  j) 

End  If 

End  If 

If  j  =  3  Then 

txtbetat(i).Text  =  beta(i,  j) 
txtXt(i).Text  =  x(i,  j) 
txtXUt(i).Text  =  Xu(i,  j) 
txtXzt(i).Text  =  Xz(i,  j) 
txtXthetat(i).Text  =  Xtheta(i,  j) 
txtXwt(i).Text  =  Xw(i,  j) 
txtrt(i).Text  =  R(i,  j) 
txtMt(i).Text  =  M(i,  j) 

If  i  =  2  Then 

txtalphat(i).Text  =  alpha(i,  j) 
txtMwt(i).Text  =  Mw(i,  j) 
txtMzt(i).Text  =  Mz(i,  j) 
txtrstt(i).Text  =  rst(i,  j) 
txtDt(i).Text  =  D(i,  j) 
txtYt(i).Text  =  Y(i,  j) 
txtYwt(i).Text  =  Yw(i,  j) 

El  self  1  =  3  Then 

txtalphat(i).Text  =  alpha(i,  j) 
txtMwt(i).Text  =  Mw(i,  j) 
txtMzt(i).Text  =  Mz(i,  j) 
txtrstt(i).Text  =  rst(i,  j) 
txtDt(i).Text  =  D(i,  j) 
txtYt(i).Text  =  Y(i,  j) 
txtYwt(i).Text  =  Yw(i,  j) 

End  If 

End  If 


txtA(i).Text  =  A(i) 
txtrht(i).Text  =  rht(i) 
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Next  j 


Next  i 

txttau.Text  =  tau 

End  Sub 

Private  Sub  UpdatePerf() 

For  i  =  1  To  3 

txtTTtl(i).Text  =  TTtl(i) 
txtPPtl(i).Text  =  PPtl(i) 

Next  i 

For  i  =  1  To  2 

txtPEPtl(i).Text  =  PEPtl(i) 

Next  i 

For  i  =  2  To  3 

txtDelta(i).Text  =  delta(i) 
txtomegasftc(i).Text  =  omegasftc(i) 
txtomegap(i).Text  =  omegap(i) 
txtomegasl(i).Text  =  omegas(i) 
txtomega(i).Text  =  omegat(i) 
txtAAl(i).Text  =  AAl(i) 
txtPhim(i).Text  =  Capphi(i,  mn) 
txtPtPtl(i).Text  =  PtPtl(i) 

Next  i 

txtTETtl(l).Text  =  TETtl(l) 
txtTRlTtl.Text  =  TrlTtl 
txtPRlPtl.Text  =  PrlPtl 
txtTt3Ttl.Text  =  TtTtl(3) 
txtPI.Text  =  PtPtl(3) 
txtB  Tau.  Text  =  tau 
txteff.Text  =  CompEff 

End  Sub 

Private  Sub  UpdateGeo() 

For  i  =  2  To  3 

txtH(i).Text  =  h(i) 
txtC(i).Text  =  C(i) 
txtZ(i).Text  =  Z(i) 
txtZrev(i).Text  =  intZ(i) 
txtARrev(i).Text  =  ARrev(i) 
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txtCrev(i).Text  =  Crev(i) 
txtS(i).Text  =  S(i) 
txtcamberh(i)  =  camber(i,  1) 
txtcamberm(i)  =  camber(i,  2) 
txtcambert(i)  =  camber(i,  3) 
txtincidenceh(i)  =  icor(i,  1) 
txtincidencem(i)  =  icor(i,  2) 
txtincidencet(i)  =  icor(i,  3) 
txtdevh(i)  =  dref(i,  1) 
txtdevm(i)  =  dref(i,  2) 
txtdevt(i)  =  dref(i,  3) 

Next  i 

End  Sub 

Private  Sub  Form_Load() 

'  Initialize  Inlet  Condition  Variables. 
alpha(l,  3)  =  fmiConipressor.txtalphat(l).Text 
gamma  =  frmCompressor.txtgamma.Text 
Mz(l,  3)  =  frmCompressor.txtMzt(l).Text 
mdot  =  frmCompressor.txtmdot.Text 
Mw(l,  3)  =  frmCompressor.txtMwt(l).Text 
omega  =  frmCompressor.txtspeed.Text 
Pt(l)  =  frmCompressor.txtPtl.Text 
Rbar  =  frmCompressor.txtR.Text 
Tt(l)  =  frmCompressor.txtTtl.Text 
'  Initialize  Rotor  Assumption  Variables 
D(2,  2)  =  frmCompressor.txtDm(2).Text 
sigma(2,  2)  =  frmCompressor.txtsigmam(2).Text 
pitch(2)  =  IrmCompressor.txtR21.Text 
phi21  =  frmCompressor.txtphi21.Text 
fsigma(2)  =  frmCompressor.txtfsigma(2).Text 
'  Initialize  Stator  Assumption  Variables 
A31  =  frmCompressor.txtA31.Text 
sigma(3,  2)  =  frmCompressor.txtsigmam(3).Text 
pitch(3)  =  frmCompressor.txtR32.Text 
phi32  =  frmCompressor.txtphi32.Text 
fsigma(3)  =  frmCompressor.txtfsigma(3).Text 


Initialize  Stage  Performance  variables 
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delta(2)  =  frmCompressor.txtDelta(2).Text 
delta(3)  =  fmiCompressor.txtDelta(3).Text 
'  initialize  blade  geometry  variables 
For  i  =  2  To  3 

For  j  =  1  To  3 

AR(i)  =  frmCompressor.txtAR(i).Text 
If  j  =  1  Then 

tc(i,  j)  =  frmCompressor.txttch(i).Text 
El  self  j  =  2  Then 

tc(i,  j)  =  frmCompressor.txttcm(i).Text 
El  self  j  =  3  Then 

tc(i,  j)  =  frmCompressor.txttct(i).Text 

End  If 

Next  j 

Next  i 

End  Sub 

Private  Sub  mnuAbout_Click() 
frmAbout.Show 

End  Sub 

Private  Sub  mnuContents_Click() 

MsgBox  'This  function  is  not  available  yet.”,  vbinformation,  "Not  Implemented" 

End  Sub 

Private  Sub  mnuExit_Click() 

End 

End  Sub 

Private  Sub  mnuOpen_Click() 
dlgPile.ShowOpen 

End  Sub 

Private  Sub  mnuPrint_Click() 
dlgPile .  ShowPrinter 
Printer.Copies  =  dlgPile. Copies 
Printer.Orientation  =  dlgPile. Orientation 
Printer.Print  A(i) 

Printer.EndDoc 

End  Sub 

Private  Sub  mnuSave_Click() 
dlgPile.  ShowSave 
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End  Sub 

Private  Sub  mnuSearch_Click() 


MsgBox  "This  function  is  not  available  yet.'',  vbinformation,  "Not  Implemented" 

End  Sub 

Public  Sub  SequenceO 
i  =  1 

InletCond 
HubCalc  i 

If  optNACA  =  True  Then 
ksh=  I 

slopegraph  =  True 

End  If 

If  optC  =  Tme  Then 
ksh=  I.I 

slopegraph  =  Ealse 

End  If 

If  optDCA  =  True  Then 
ksh  =  0.7 

slopegraph  =  Ealse 

End  If 

Eor  i  =  2  To  3 

MeanCalc  i 

Next  i 

loopcount  =  I 
Do 

Eor  i  =  2  To  3 

If  loopcount  =  1  Then 
omegat(i)  =  0 
omegasftc(i)  =  0 
omegas  (i)  =  0 

End  If 

tempvar(i)  =  omegat(i) 

If  i  =  2  Then 

RotorPerf  i,  loopcount 
El  self  i  =  3  Then 

StatorPerf  i,  loopcount 

End  If 
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HubTipCalc  i 
BladeGeo  i,  loopcount 


Next  i 

loopcount  =  loopcount  +  1 

Loop  Until  Abs(omegat(2)  -  tempvar(2))  <  testvar  And  Abs(omegat(3)  -  tempvar(3))  < 
testvar 

For  i  =  2  To  3 

For  j  =  1  To  3 

Incidence  i,  j,  slopegraph,  ksh 

Next  j 

Next  i 

CompEff  =  ((PtPtl(3)  *  glg(gamma))  -  1)  /  (tau  -  1) 

End  Sub 

Private  Sub  txtA31_KeyPress(Key Ascii  As  Integer) 

Select  Case  Key  Ascii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtA31.Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtalphat_KeyPress(Index  As  Integer,  KeyAscii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtalphat(l).Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 
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End  If 


Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtAR_KeyPress(Index  As  Integer,  KeyAscii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtAR(2).Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

Elself  Len(txtAR(3).Text)  <>  0  Then 
KeyAscii  =  0 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtDelta_KeyPress(Index  As  Integer,  KeyAscii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtDelta(2).Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

Elself  Len(txtDelta(3).Text)  <>  0  Then 
KeyAscii  =  0 
Beep 

End  If 
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Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtDm_KeyPress(Index  As  Integer,  KeyAscii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtDm(2).Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtgamma_KeyPress(Key Ascii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtgamma.Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtmdot_KeyPress(Key Ascii  As  Integer) 

Select  Case  KeyAscii 
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Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtmdot.Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtMwt_KeyPress(Index  As  Integer,  KeyAscii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtMwt(I).Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtMzt_KeyPress(Index  As  Integer,  KeyAscii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtMzt(l).Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 
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End  If 


Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtphi21_KeyPress(Key Ascii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtphi2I.Text) <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtphi32_KeyPress(Key Ascii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtphi32.Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtPtl_KeyPress(Key Ascii  As  Integer) 
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Select  Case  Key  Ascii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtPtl.Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtR_KeyPress(KeyAscii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtR.Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtR21_KeyPress(Key Ascii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtR21.Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
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Beep 


End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtR32_KeyPress(KeyAscii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtR32.Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtsigmam_KeyPress(Index  As  Integer,  KeyAscii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtsigmam(2).Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

Elself  Len(txtsigmam(3).Text)  <>  0  Then 
KeyAscii  =  0 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
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Beep 

End  Select 

End  Sub 

Private  Sub  txtspeed_KeyPress(Key Ascii  As  Integer) 

Select  Case  Key  Ascii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtspeed.Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txttch_KeyPress(Index  As  Integer,  KeyAscii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txttch(2).Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

Elself  Len(txttch(3).Text)  <>  0  Then 
KeyAscii  =  0 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txttcm_KeyPress(Index  As  Integer,  KeyAscii  As  Integer) 
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Select  Case  Key  Ascii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txttcm(2).Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

Elself  Len(txttcm(3).Text)  <>  0  Then 
KeyAscii  =  0 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txttct_KeyPress(Index  As  Integer,  KeyAscii  As  Integer) 

Select  Case  KeyAscii 

Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txttct(2).Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

Elself  Len(txttct(3). Text)  <>  0  Then 
KeyAscii  =  0 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

Private  Sub  txtTtl_KeyPress(KeyAscii  As  Integer) 

Select  Case  KeyAscii 
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Case  vbKeyO  To  vbKey9 
Case  vbDecimal,  46 

Case  vbKeyBack,  vbKeyClear,  vbKeyTab,  vbKeyUp,  vbKeyDown 
Case  45 

If  Len(txtTtl.Text)  <>  0  Then 

KeyAscii  =  0  '  ignore  keystroke 
Beep 

End  If 
Case  Else 

KeyAscii  =  0 
Beep 

End  Select 

End  Sub 

5.  Blade  Number  Input  Screen  Code 

Option  Explicit 
Option  Base  1 
Private  x  As  Variant 
Private  Sub  OKButton_Click() 

X  =  txtZ.Text 

If  IsNumeric(x)  =  Ealse  Then 

MsgBox  "Input  is  not  a  number!",  vbExclamation,  "Numeric  Validation" 

Else 

frmZ.Hide 

End  If 

End  Sub 

6.  About  Screen  Code 

Option  Explicit 
'  Reg  Key  Security  Options... 

Const  READ_CONTROL  =  &H20000 

Const  KEY_QUERY_VALUE  =  &H1 

Const  KEY_SET_VALUE  =  &H2 

Const  KEY_CREATE_SUB_KEY  =  &H4 

Const  KEY_ENUMERATE_SUB_KEYS  =  &H8 

Const  KEY_NOTIEY  =  &H10 

Const  KEY_CREATE_LINK  =  &H20 
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Const  KEY_ALL_  ACCESS  =  KEY_QUERY_VALUE  +  KEY_SET_VALUE  + 

_  KEY_CREATE_SUB_KEY  +  KEY_ENUMERATE_SUB_KEYS  + 

_  KEY_NOTIEY  +  KEY_CREATE_LINK  +  READ_CONTROL 
'  Reg  Key  ROOT  Types... 

Const  HKEY_LOCAL_MACHINE  =  &H80000002 
Const  ERROR_SUCCESS  =  0 

Const  REG_SZ  =  1  '  Unicode  nul  terminated  string 

Const  REG_DWORD  =  4  '  32-bit  number 

Const  gREGKEYSYSINEOLOC  =  "SOFTWARE\Microsoft\Shared  Tools  Location" 

Const  gREGVALSYSINFOLOC  =  "MSINFO” 

Const  gREGKEYSYSINFO  =  "SOFTWARE\Microsoft\Shared  Tools  \MSINFO'' 

Const  gREGVALSYSINFO  =  "PATH” 

Private  Declare  Function  RegOpenKeyEx  Lib  ”advapi32"  Alias  "RegOpenKeyExA"  (ByVal  hKey 
As  Long,  ByVal  IpSubKey  As  String,  ByVal  ulOptions  As  Long,  ByVal  samDesired  As 
Long,  ByRef  phkResult  As  Long)  As  Long 

Private  Declare  Function  RegQueryValueEx  Lib  ”advapi32”  Alias  "RegQueryValueExA"  (ByVal 
hKey  As  Long,  ByVal  IpValueName  As  String,  ByVal  IpReserved  As  Long,  ByRef 
IpType  As  Long,  ByVal  IpData  As  String,  ByRef  IpcbData  As  Long)  As  Long 

Private  Declare  Function  RegCloseKey  Lib  ”advapi32"  (ByVal  hKey  As  Long)  As  Long 

Private  Sub  cmdSysInfo_Click() 

Call  StartSysInfo 

End  Sub 

Private  Sub  cmdOK_Click() 

Unload  Me 

End  Sub 

Private  Sub  Form_Load() 

Me. Caption  =  "About  ”  &  App. Title 

IblVersion. Caption  =  "Version  ”  &  App. Major  &  &  App. Minor  &  &  App. Revision 

IblTitle. Caption  =  App. Title 

IblDescription. Caption  =  App.FileDescription 

End  Sub 

Public  Sub  StartSysInfoO 

On  Error  GoTo  SysInfoErr 

Dim  rc  As  Long 

Dim  SysInfoPath  As  String 

'  Try  To  Get  System  Info  Program  Path\Name  From  Registry... 
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gREGKEYSYSINFO, 


If  GetKeyValue(HKEY_LOCAL_MACHINE, 

gREGVALSYSINEO,  SysInfoPath)  Then 

'  Try  To  Get  System  Info  Program  Path  Onfy  From  Registry... 

Elself  GetKeyValue(HKEY_LOCAL_MACHINE,  gREGKEYSYSINFOLOC, 

gREGVALSYSINFOLOC,  SysInfoPath)  Then 

'  Validate  Existance  Of  Known  32  Bit  File  Version 

If  (Dir(SysInfoPath  &  ''\MSINF032.EXE'')  <>  Then 

SysInfoPath  =  SysInfoPath  &  ’'\MSINF032.EXE' 

'  Error  -  File  Can  Not  Be  Found... 

Else 

GoTo  SysInfoErr 

End  If 

'  Error  -  Registry  Entry  Can  Not  Be  Found... 

Else 

GoTo  SysInfoErr 

End  If 

Call  Shell(SysInfoPath,  vbNormalFocus) 

Exit  Sub 
SysInfoErr: 

MsgBox  "System  Information  Is  Unavailable  At  This  Time”,  vbOKOnly 

End  Sub 

Public  Function  GetKeyValue(KeyRoot  As  Long,  KeyName  As  String,  SubKeyRef  As  String, 
ByRef  KeyVal  As  String)  As  Boolean 

Dim  i  As  Long  '  Loop  Counter 

Dim  rc  As  Long  '  Return  Code 

Dim  hKey  As  Long  '  Handle  To  An  Open  Registry  Key 

Dim  hDepth  As  Long 

Dim  KeyValType  As  Long  '  Data  Type  Of  A  Registry  Key 

Dim  tmpVal  As  String  '  Tempory  Storage  For  A  Registry  Key  Value 

Dim  KeyValSize  As  Long  '  Size  Of  Registry  Key  Variable 


'  Open  RegKey  Under  KeyRoot  {HKEY_LOCAL_MACHINE... } 


rc  =  RegOpenKeyEx(KeyRoot,  KeyName,  0,  KEY_ALL_ACCESS,  hKey) 
'  Open  Registry  Key 

If  (rc  <>  ERROR_SUCCESS)  Then  GoTo  GetKeyError 
tmpVal  =  String$(I024,  0) 

KeyValSize  =  1024 
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Handle  Error... 

Allocate  Variable  Space 
Mark  Variable  Size 


'  Retrieve  Registry  Key  Value... 


rc  =  RegQueryValueEx(hKey,  SubKeyRef,  0,  _ 

KeyValType,  tmpVal,  KeyValSize)  '  Get/Create  Key  Value 

If  (rc  <>  ERROR_SUCCESS)  Then  GoTo  GetKeyError  '  Handle  Errors 

If  (Asc(Mid(tmpVal,  KeyValSize,  1))  =  0)  Then  '  Win95  Adds  Null 

Terminated  String... 

tmpVal  =  Left(tmpVal,  KeyValSize  -  1)  '  Null  Eound,  Extract 

Erom  String 

Else  '  WinNT  Does  NOT  Null  Terminate  String... 

tmpVal  =  Left(tmpVal,  KeyValSize)  '  Null  Not  Eound, 

Extract  String  Only 

End  If 


'  Determine  Key  Value  Type  Eor  Conversion... 


Select  Case  KeyValType 
Case  REG_SZ 

KeyVal  =  tmpVal 
Case  REG_DWORD 

Eor  i  =  Len(tmpVal)  To  1  Step-1 


Search  Data  Types... 

String  Registry  Key  Data  Type 
Copy  String  Value 

Double  Word  Registry  Key  Data  Type 
Convert  Each  Bit 

Build  Value  Char. 


KeyVal  =  KeyVal  +  Hex(Asc(Mid(tmpVal,  i,  1))) 
By  Char. 


Next 


KeyVal  =  Eormat$("&h"  +  KeyVal) 
End  Select 
GetKey  Value  =  True 
rc  =  RegCloseKey(hKey) 

Exit  Eunction 
GetKeyError: 

KeyVal  = 

GetKey  Value  =  Ealse 
rc  =  RegCloseKey(hKey) 

End  Eunction 


'  Convert  Double  Word  To  String 


'  Return  Success 
'  Close  Registry  Key 
'  Exit 

Cleanup  After  An  Error  Has  Occured... 
Set  Return  Val  To  Empty  String 
Return  Eailure 
Close  Registry  Key 


7. 


Module  mdllnletCond  Code 


Option  Explicit 
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Option  Base  1 

Public  gamma  As  Double 

'  Specific  Heat  Ratio  (Cp/Cv) 

Public  mdot  As  Double 
'  Mass  Flow 
Public  Rbar  As  Double 
'  Specific  Gas  Constant 
Public  omega  As  Double 
'  wheel  speed 
Public  Cp  As  Double 
'  at  constant  pressure 
Const  g  As  Long  =  32.2 
'  gravitational  constant 
Public  Capphi(3,  3)  As  Double 
'  Flow  Function 
Public  Tt(3)  As  Double 
'  Total  Temperature 
Public  Pt(3)  As  Double 
'  Total  Pressure 
Public  M(3,  3)  As  Double 
'  Mach  Number 
Public  Mw(3,  3)  As  Double 
'  Mach  Relative  to  the  Blade 
Public  Mz(3,  3)  As  Double 
'  Mach  of  the  Axial  Component 
Public  alpha(3,  3)  As  Double 
'  Inlet  Flow  angle 
Public  beta(3,  3)  As  Double 
'  Inlet  Flow  angle 
Public  x(3,  3)  As  Double 
'  Dimensionless  Velocity 
Public  Xz(3,  3)  As  Double 

'  Dimensionless  Velocity  of  the  Axial  Component 
Public  Xtheta(3,  3)  As  Double 
'  Dimensionaless  Velocity  along  Theta 
Public  Xu(3,  3)  As  Double 

'  Dimensionless  Velocity  of  the  Wheel  Speed  Component 
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Public  Xw(3,  3)  As  Double 
'  Dimensionless  Velocity  Relative  to  the  Blade 
Public  rot(3)  As  Double 
'  Density 

Public  Vt(3)  As  Double 
'  Total  Velocity 
Public  A(3)  As  Double 
'  Annulus  Area 
Public  R(3,  3)  As  Double 
'  Radius 

Public  rht(3)  As  Double 
'  Hub  to  Tip  Ratio 
Public  pitch(3)  As  Double 
'  Mean  Line  Pitch  (Rm2/Rml) 

Public  rst(3,  3)  As  Double 

'  Degree  of  Reaction 

Public  i  As  Integer 

Public  Const  hub  As  Integer  =  I 

'  Integer  used  when  making  hub  calculations  (hub  =1) 

Public  Const  mn  As  Integer  =  2 

'  Integer  used  when  making  mean  line  calculations  (mn  =  2) 

Public  Const  tip  As  Integer  =  3 
'  Integer  used  when  making  tip  calculations  (tip  =  3) 

Public  Sub  InletCondO 

'  Calculation  of  Inlet  Conditions  at  the  tip 
i  =  1 

beta(i,  tip)  =  Arccos(Mz(i,  tip)  /  Mw(i,  tip)) 

M(i,  tip)  =  Mz(i,  tip)  /  Cos(DegToRad(alpha(i,  tip))) 

x(i,  tip)  =  Sqr((((gamma  -  1)  /  2)  *  M(i,  tip)  *  2)  /  (1  +  (((gamma  -  1)  /  2)  *  M(i,  tip)  2))) 
Xz(i,  tip)  =  x(i,  tip)  *  Cos(DegToRad(alpha(i,  tip))) 

Xtheta(i,  tip)  =  Xz(i,  tip)  *  Tan(DegToRad(alpha(i,  tip))) 

Xu(i,  tip)  =  Xtheta(i,  tip)  -i-  Xz(i,  tip)  *  Tan(DegToRad(beta(i,  tip))) 

Xw(i,  tip)  =  Xwfunc(Xz(i,  tip),  beta(i,  tip)) 
rot(i)  =  Pt(i)  /  (12  *  Rbar  *  Tt(i)) 

Cp  =  gg  1  (gamma)  *  Rbar 

Vt(i)  =  Sqr(2  *  Cp  *  g  *  Tt(i))  *  12 

Capphi(i,  tip)  =  x(i,  tip)  *  (1  -  x(i,  tip)  *  2)  *  (1  /  (gamma  -  1)) 
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A(i)  =  (mdot  /  (rot(i)  *  Vt(i)))  *  (1  /  (Capphi(i,  tip)  *  Cos(DegToRad(alpha(i,  tip))))) 

R(i,  tip)  =  (Xu(i,  tip)  *  Vt(i))  /  omega 

R(i,  hub)  =  Sqr(R(i,  tip)  2  -  A(i)  /  (22  /  7)) 

rht(i)  =  rhtfunc(R(i,  tip),  R(i,  hub)) 

'  Mean  Calculations 

R(i,  mn)  =  rmfunc(R(i,  tip),  R(i,  hub)) 

Xtheta(i,  mn)  =  Xthetafunc(R(i,  tip),  R(i,  mn),  Xtheta(i,  tip)) 

Xz(i,  mn)  =  Xz(i,  tip) 

alpha(i,  mn)  =  alphafunc(Xtheta(i,  mn),  Xz(i,  mn)) 

Xu(i,  mn)  =  Xufunc(R(i,  tip),  R(i,  mn),  Xu(i,  tip)) 
beta(i,  mn)  =  betafunc(Xu(i,  mn),  Xtheta(i,  mn),  Xz(i,  mn)) 
x(i,  mn)  =  Xfunc(Xz(i,  mn),  alpha(i,  mn)) 

M(i,  mn)  =  Mach(gamma,  x(i,  mn)) 

Xw(i,  mn)  =  Xwfunc(Xz(i,  mn),  beta(i,  mn)) 

Mw(i,  mn)  =  Mach(gamma,  Xw(i,  mn)) 

Mz(i,  mn)  =  Machz(M(i,  mn ),  alpha(i,  mn)) 

End  Sub 

8.  Module  mdlHubCalc  Code 

Option  Explicit 
Option  Base  1 

Public  sigma(3,  3)  As  Double 
'  Solidity 

Public  Sub  HubCalc(i  As  Integer) 

'  Hub  Calculations 

Xz(i,  hub)  =  Xz(i,  mn)  '  Radial  Equilibrium 

Xtheta(i,  hub)  =  Xthetafunc(R(i,  mn),  R(i,  hub),  Xtheta(i,  mn)) 

alpha(i,  hub)  =  alphafunc(Xtheta(i,  hub),  Xz(i,  hub)) 

x(i,  hub)  =  Xfunc(Xz(i,  hub),  alpha(i,  hub)) 

Xu(i,  hub)  =  Xufunc(R(i,  mn),  R(i,  hub),  Xu(i,  mn)) '  for  constant  wheel  speed 
beta(i,  hub)  =  betafunc(Xu(i,  hub),  Xtheta(i,  hub),  Xz(i,  hub)) 

Xw(i,  hub)  =  Xwfunc(Xz(i,  hub),  beta(i,  hub)) 

M(i,  hub)  =  Mach(gamma,  x(i,  hub)) 

Mw(i,  hub)  =  Mach(gamma,  Xw(i,  hub)) 

Mz(i,  hub)  =  Machz(M(i,  hub),  alpha(i,  hub)) 

End  Sub 
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9. 


Module  mdlMeanCalc  Code 


Option  Explicit 
Option  Base  1 
Public  Y(3,  3)  As  Double 
'  Dimensionless  Local  Velocity 
Public  Yw(3,  3)  As  Double 

'  Dimensionless  Local  Velocity  Relative  to  the  Blade 
Public  Yz(3,  3)  As  Double 

'  Dimensionless  Local  Velocity  along  the  Axial  Component 
Public  A3 1  As  Double 
'  Ratio  of  alpha.  Unity  for  a  repeating  stage. 

Public  D(3,  3)  As  Double 

'  Diffusion  Factor 

Public  tau  As  Double 

'  Total  temperature  ratio 

Public  phi21  As  Double 

'  Change  in  Axial  Velocity  across  the  Rotor. 

Public  phi32  As  Double 
'  Change  in  Axial  Velocity  across  the  Stator. 

Public  phi2m  As  Double 

'  Change  in  Axial  Velocity  at  the  Rotor  exit  at  the  Mean  Line. 

Public  philm  As  Double 
'  ratio  of  axial  velocity  to  rotation  velocity 
Public  Sub  MeanCalc(i  As  Integer) 

'  Mean  Calculations  for  the  Rotor  or  Stator 
If  i  =  2  Then 

Dim  temp  As  Double 

beta(i,  rm)  =  Arcsin(SinB2(DiffA(sigma(i,  mn),  pitch(i),  D(i,  mn),  phi21,  beta(i 
- 1,  mn)),  DiffB(sigma(i,  mn),  pitch(i)),  pitch(i))) 

Xu(i,  mn)  =  pitch(i)  *  Xu(i  -  1,  mn) 

philm  =  Xz(i  -  1,  mn)  /  Xu(i  -  1,  mn) 

phi2m  =  phi21  *  philm  *  (1  /  pitch(i)) 

Xz(i,  mn)  =  phi2m  *  Xu(i,  mn) 

Xtheta(i,  mn)  =  Xu(i,  mn)  -  Xz(i,  mn)  *  Tan(DegToRad(beta(i,  mn))) 

alpha(i,  mn)  =  alphafunc(Xtheta(i,  mn),  Xz(i,  mn)) 

tau  =  taufunc(Xu(i,  mn),  Xtheta(i,  mn),  Xu(i  -  1,  mn),  Xtheta(i  -  1,  mn)) 

Else 
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alpha(i,  mn)  =  A3 1  *  alpha(i  -  2,  mn) 

D(i,  mn)  =  Diffusion(phi32,  alpha(i  -  1,  mn),  alpha(i,  mn),  pitch(i),  sigma(i, 
mn)) 

Xz(i,  mn)  =  phi32  *  Xz(i  -  1,  mn) 

Xtheta(i,  mn)  =  Xz(i,  mn)  *  Tan(DegToRad(alpha(i,  mn))) 

Xu(i,  mn)  =  Xu(i  -  1,  mn)  *  pitch(i) 

beta(i,  mn)  =  betafunc(Xu(i,  mn),  Xtheta(i,  mn),  Xz(i,  mn)) 

End  If 

rst(i,  mn)  =  DofReaction(Xtheta(i  -  1,  mn),  Xu(i  -  1,  mn),  Xtheta(i,  mn),  Xu(i,  mn)) 
x(i,  mn)  =  Xfunc(Xz(i,  mn),  alpha(i,  mn)) 

Xw(i,  mn)  =  Xwfunc(Xz(i,  mn),  beta(i,  mn)) 

Y(i,  mn)  =  Yfunc(x(i,  mn),  tau) 

Yw(i,  mn)  =  Yfunc(Xw(i,  mn),  tau) 

M(i,  mn)  =  Mach(gamma,  Y(i,  mn)) 

Mw(i,  mn)  =  Mach(gamma,  Yw(i,  mn)) 

Mz(i,  mn)  =  Machz(M(i,  mn),  alpha(i,  mn)) 

End  Sub 

10.  Module  mdlStgPerformance 

Option  Explicit 
Option  Base  1 

Public  omegasftc(3)  As  Double 
'  Secondary  flow  and  tip  clearance  loss 
Public  omegap(3)  As  Double 
'  Profile  Loss 

Public  omegas(3)  As  Double 
'  Shock  Loss 

Public  omegat(3)  As  Double 

'  Total  Loss  (secondary  flow  +  tip  clearance  +  profile) 

Public  TTtl(3)  As  Double 
'  Static  to  Total  Temperature  Ratio 
Public  PPtl(3)  As  Double 
'  Static  to  Total  Pressure  Ratio 
Public  TrlTtl  As  Double 
'  Total  Relative  Temperature  Ratio 
Public  PrlPtl  As  Double 
'  Total  Relative  Pressure  Ratio 
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Public  TETtl(3)  As  Double 
Public  PEPtl(3)  As  Double 
Public  PtPtl(3)  As  Double 
'  Total  Pressure  Ratio 
Public  TtTtl(3)  As  Double 
'  Total  Temperature  Ratio 
Public  AA1(3)  As  Double 
'  Area  constriction  ratio 
Public  delta(3)  As  Double 
'  Tip  Gap 

Public  Sub  RotorPerf(i  As  Integer,  loopcount  As  Integer) 

Dim  ShLoss  As  Double 

omegap(i)  =  2  *  sigma(i,  mn)  *  ((Cos(DegToRad(beta(i  -  1,  mn)))  ^  1)  ! 
(Cos(DegToRad(beta(i,  mn)))  3))  *  (0.005  +  (0.16  *  (D(i,  mn)  4))) 

If  loopcount  >  1  Then 

If  Mw(i  -  1,  mn)  >  1  Then 

ShLoss  =  ((1  -  ShockLoss(Mw(i  -  1,  mn),  gamma))  /(!-(!+  ((gamma 
- 1)  /  2)  *  Mw(i-  1,  mn)  2)  (-1  *  ggl(gamma)))) 

End  If 

omegas(i)  =  ShLoss 

omegasftc(i)  =  SFTC(beta(i  -  1,  mn),  beta(i,  mn),  sigma(i,  mn),  h(i),  delta(i), 
S(i)) 

Else 

omegas  (i)  =  0 
omegasftc(i)  =  0 

End  If 

omegat(i)  =  omegap(i)  +  omegasftc(i)  +  omegas(i) 

If  loopcount  =  1  Then 

TTtl(i-  1)=  1  -(x(i-  l,mn)''2) 

PPtl(i  -  1)  =  TTtl(i  -  1)  ggl(gamma) 

TrlTtl  =  TTtl(i  -  1)  +  (Xw(i  -  1,  mn)  2) 

PrlPtl  =  TrlTtl  *  ggl(gamma) 

TETtl(i  -  1)  =  TrlTtl  +  (Xu(i,  mn)  2)  -  (Xu(i  -  1,  mn)  2) 

PEPtl(i  -  1)  =  TETtl(i  -  1)  ggl(gamma) 

End  If 

PEPtl(i)  =  PEPtl(i  -  1)  -  (omegat(i)  *  (PrlPtl  -  PPtl(i  -  1))) 

TtTtl(i)  =  tau 
TETtl(i)  =  TETtl(i-  1) 
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PtPtl(i)  =  PEPtl(i)  *  ((TtTtl(i)  /  TETtl(i))  *  ggl(gamma)) '  Rotor  Compression  ratio 
TTtl(i)  =  TtTtl(i)  -  (x(i,  mn)  2) 

PPtl(i)  =  PtPtl(i)  *  ((TTtl(i)  /  tau)  ggl(gamma)) 

Capphi(i,  mn)  =  Y(i,  mn)  *  ((1  -  (Y(i,  mn)  2))  gl  (gamma)) 

AAl(i)  =  ((Capphi(i  -  1,  tip)  *  Cos(DegToRad(alpha(i  -  1,  tip))))  /  (Capphi(i,  mn)  * 
Cos(DegToRad(alpha(i,  mn)))))  *  ((Sqr(TtTtl(i)))  /  PtPtl(i)) 

A(i)  =  AAl(i)  *  A(i  -  1) 

R(i,  mn)  =  R(i  -  1,  mn)  *  pitch(i) 

h(i)  =  A(i)  /  (2  *  (22  /  7)  *  R(i,  mn)) 

rht(i)  =  rhtfunc2(h(i),  R(i,  mn)) 

R(i,  tip)  =  (2  /  (1  +  rht(i)))  *  R(i,  mn) 

R(i,  hub)  =  rht(i)  *  R(i,  tip) 

End  Sub 


Public  Sub  StatorPerf(i  As  Integer,  loopcount  As  Integer) 

Dim  ShLoss  As  Double 

omegap(i)  =  2  *  sigma(i,  mn)  *  ((Cos(DegToRad(alpha(i  -  1,  mn)))  *  2)  / 
(Cos(DegToRad(alpha(i,  mn)))  ^  3))  *  (0.005  +  (0.16  *  (D(i,  mn)  4))) 

If  loopcount  >  1  Then 

If  M(i  -  I ,  mn)  >  I  Then 

ShLoss  =  ((I  -  ShockLoss(M(i  -  I,  mn),  gamma))  /  (I  -  (I  +  ((gamma  - 
I)  /  2)  *  M(i  -  I,  mn)  *  2)  *  (-1  *  ggl(gamma)))) 

End  If 

omegas(i)  =  ShLoss 

omegasftc(i)  =  SFTC(alpha(i  -  I,  mn),  alpha(i,  mn),  sigma(i,  mn),  h(i),  delta(i), 
S(i)) 

Else 

omegas  (i)  =  0 
omegasftc(i)  =  0 

End  If 

omegat(i)  =  omegap(i)  +  omegasftc(i)  +  omegas(i) 

PtPtl(i)  =  PtPtI(i  -  I)  -  omegat(i)  *  (PtPtI(i  -  I)  -  PPtI(i  -  I)) 

TtTtl(i)  =  tau 

TTtl(i)  =  TtTtl(i)  -  (x(i,  mn)  2) 

PPtl(i)  =  ((TTtl(i)  /  TtTtl(i))  ggl(gamma))  *  PtPtl(i) 

Capphi(i,  mn)  =  Y(i,  mn)  *  ((I  -  (Y(i,  mn)  2))  *  g  I  (gamma)) 

AAI(i)  =  ((Capphi(i  -  2,  tip)  *  Cos(DegToRad(alpha(i  -  2,  tip))))  /  (Capphi(i,  mn)  * 
Cos(DegToRad(alpha(i,  mn)))))  *  ((Sqr(TtTtI(i)))  /  PtPtl(i)) 


A(i)  =  AAI(i)*A(i-2) 
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R(i,  mn)  =  R(i  -  1,  mn)  *  pitch(i) 
h(i)  =  A(i)  /  (2  *  (22  /  7)  *  R(i,  mn)) 
rht(i)  =  rhtfunc2(h(i),  R(i,  mn)) 

R(i,  tip)  =  (2  /  (1  +  rht(i)))  *  R(i,  mn) 

R(i,  hub)  =  rht(i)  *  R(i,  tip) 

End  Sub 

11.  Module  mdlHubTipCalc  Code 

Option  Explicit 
Option  Base  1 
Dim  1  As  Integer 
Public  fsigma(3)  As  Double 

'  Linear  function  of  sigma 

Public  Sub  HubTipCalc(i  As  Integer) 

'  Hub  &  Tip  Calculations  for  the  Rotor  or  Stator 
Eor  1  =  1  To  3  Step  2 

Xz(i,  1)  =  Xz(i,  mn) 

Xtheta(i,  1)  =  (R(i,  mn)  *  Xtheta(i,  mn))  /  R(i,  1) '  Radial  Equilibrium 
Xu(i,  1)  =  Xufunc(R(i,  mn),  R(i,  1),  Xu(i,  mn))  'for  constant  wheel  speed 
alpha(i,  1)  =  alphafunc(Xtheta(i,  1),  Xz(i,  1)) 
beta(i,  1)  =  betafunc(Xu(i,  1),  Xtheta(i,  1),  Xz(i,  1)) 
x(i,  1)  =  Xfunc(Xz(i,  1),  alpha(i,  1)) 

Xw(i,  1)  =  Xwfunc(Xz(i,  1),  beta(i,  1)) 

Y(i,  1)  =  Yfunc(x(i,  1),  tau) 

Yw(i,  1)  =  Yfunc(Xw(i,  1),  tau) 

M(i,  1)  =  Mach(gamma,  Y(i,  1)) 

Mw(i,  1)  =  Mach(gamma,  Yw(i,  1)) 

Mz(i,  1)  =  Machz(M(i,  1),  alpha(i,  1)) 

If  1  =  1  Then 

sigma(i,  1)  =  sigma(i,  mn)  *  (2  -  fsigma(i)) 

El  self  1  =  3  Then 

sigma(i,  1)  =  sigma(i,  mn)  *  fsigma(i) 

End  If 
If  i  =  2  Then 

D(i,  1)  =  Diffusion(phi2I,  beta(i  -  1,1),  beta(i,  1),  pitch(i),  sigma(i,  1)) 

Else 


64 


D(i,  1)  =  Diffusion(phi32,  alpha(i-  1, 1),  alpha(i,  1),  pitch(i),  sigma(i,  1)) 

End  If 

rst(i,  1)  =  DofReaction(Xtheta(i  -  1,1),  Xu(i  -  1,1),  Xtheta(i,  1),  Xu(i,  1)) 

Next  1 

End  Sub 

12.  Module  mdlBladeGeometry  Code 

Option  Explicit 
Option  Base  1 
Public  AR(3)  As  Double 
'  Aspect  Ratio 
Public  C(3)  As  Double 
'  Chord 

Public  Crev(3)  As  Double 
'  revised  chord  length 
Public  ARrev(3)  As  Double 
'  revised  aspect  ratio 
Public  Z(3)  As  Double 
'  Number  of  Blades  calculated 
Public  intZ(3)  As  Integer 
'  Number  of  Blades  chosen  by  the  user. 

Public  h(3)  As  Double 
'  Blade  height 
Public  S(3)  As  Double 
'  blade  spacing 

Public  Sub  BladeGeo(i  As  Integer,  loopcount  As  Integer) 
h(i)  =  R(i,  tip)  -  R(i,  hub) '  blade  height 
C(i)  =  h(i)  /  AR(i) '  Calculated  chord  length 

Z(i)  =  2  *  (22  /  7)  *  R(i,  mn)  *  sigma(i,  mn)  /  C(i) '  Calculated  no.  of  blades 
If  loopcount  =  1  Then 
If  i  =  2  Then 

frmZ.lblZ. Caption  =  "The  calculated  no.  of  blades  in  the  rotor  (Zr)  is 
&Z(i) 

El  self  1  =  3  Then 

frmZ.lblZ. Caption  =  "The  calculated  no.  of  blades  in  the  stator  (Zs)  is 
&Z(i) 

End  If 

frmZ.Show  vbModal 
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intZ(i)  =  frmZ.txtZ.Text 


End  If 

ARrev(i)  =  h(i)  /  ((2  *  (22  /  7)  *  R(i,  mn)  *  sigma(i,  mn))  /  intZ(i)) '  Revised  aspect  ratio 
frmZ.txtZ.Text  = 

Crev(i)  =  h(i)  /  ARrev(i) '  Revised  chord  length 
S(i)  =  Crev(i)  /  sig  ma(i,  mn) '  Blade  spacing 

End  Sub 

13.  Module  mdlincidence  Code 

Option  Explicit 
Option  Base  1 

Public  i0(6,  9)  As  Double,  slopen(7,  9)  As  Double,  ikit(6)  As  Double,  i0delta0(6,  9)  As  Double, 
slopemn(7,  9)  As  Double,  Kit(7,  7)  As  Double,  d0(6,  9)  As  Double,  slopem(7,  9)  As 
Double,  dkit(6)  As  Double 

'  X  values  for  6th  order  curve  fit  polynomial. 

Public  i010(3,  3)  As  Double 

'  Zero -camber  incidence  angle 

Public  d010(3,  3)  As  Double 

'  Zero -camber  deviation  angle 

Public  N(3,  3)  As  Double 

'  incidence  angle  slope  factor  n 

Public  ikt(3,  3)  As  Double 

'  incidence  angle  correction  factor  for  thickness 

Public  dkt(3,  3)  As  Double 

'  deviation  angle  correction  factor  for  thickness 

Public  i0ref(3,  3)  As  Double 

'  reference  incidence  angle 

Public  d0ref(3,  3)  As  Double 

'  reference  deviation  angle 

Public  i2d  As  Double 

'  2  dimentional  incidence  angle 

Public  icor(3,  3)  As  Double 

'  corrected  2  dimentional  incidence  angle 

Public  dref(3,  3)  As  Double 

'  deviation  angle 

Public  camber(3,  3)  As  Double 

'  camber  angle 

Public  i0d0(3,  3)  As  Double 
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'  variation  of  iO  -  dO 

Public  ktbar(3,  3)  As  Double 

'  camber  angle  correction  factor  for  thickness 

Public  onemn(3,  3)  As  Double 

'  camber  angle  slope  factor  1  -m+n 

Public  dm(3,  3)  As  Double 

'  deviation  angle  slope  factor  m 

Public  tc(3,  3)  As  Double 

'  thickness  to  cord  ratio. 

Public  k  As  Integer 
Dim  o  As  Integer,  p  As  Single 
Public  Xvar()  As  Double 
Public  Yvar()  As  Double 

Public  Sub  Incidence(i  As  Integer,]  As  Integer,  slopegraph  As  Boolean,  ksh  As  Single) 
'  Variation  of  (i0)10-(delta0)10 
'  Load  constants 
'  solidity  =  0.4 

i0delta0(6,  1)  = -1.72244891587855E-1 1 
i0delta0(5,  1)  =  3.05183134666209E-09 
i0delta0(4,  1)  =  -2.54582401992831E-07 
iOdeltaO(3,  1)  =  8.9308266844057E-06 
i0delta0(2,  1)  = -1.64583732868095E-04 
i0delta0(l,  1)  =  2.36162560875073E-02 
'  solidity  =  0.6 

i0delta0(6,  2)  =  -1.21461996934098E-10 
i0delta0(5,  2)  =  2.44632274559037E-08 
i0delta0(4,  2)  =  -1.89083802948353E-06 
i0delta0(3,  2)  =  6.5 102496755287 lE-05 
i0delta0(2,  2)  =  -1.10036026944726E-03 
i0delta0(l,  2)  =  4.37180120798075E-02 
'  solidity  =  0.8 

i0delta0(6,  3)  =  -1.56214579635869E-10 
i0delta0(5,  3)  =  3.09609686430234E-08 
i0delta0(4,  3)  =  -2.3589555206982E-06 
i0delta0(3,  3)  =  7.920583 15148836E-05 
i0delta0(2,  3)  =  -1.27482005365209E-03 
i0delta0(l,  3)  =  5.69282884198401E-02 
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'  solidity  =1.0 

i0delta0(6,  4)  =  -6.93736034940097E-1 1 
i0delta0(5,  4)  =  8.79693724809005E-09 
i0delta0(4,  4)  =  -3.60272149557694E-07 
i0delta0(3,  4)  =  -4.42503626629787E-07 
i0delta0(2,  4)  =  1.21654421718631E-05 
i0delta0(l,  4)  =  6.527324366 19772E-02 
'  solidity  =1.2 

i0delta0(6,  5)  =  1.M197027145803E-10 
i0delta0(5,  5)  =  -2.10810745440021E-08 
i0delta0(4,  5)  =  1. 37089527 130208E-06 
i0delta0(3,  5)  =  -3.80616685049517E-05 
i0delta0(2,  5)  =  8.3487583 1969839E-05 
i0delta0(l,  5)  =  8.43313899049463E-02 
'  solidity  =1.4 

i0delta0(6,  6)  =  5.17847890355591E-11 
i0delta0(5,  6)  =  -1.11055450426056E-08 
i0delta0(4,  6)  =  6.341 13 19506379E-07 
i0delta0(3,  6)  =  -1.51366527276764E-05 
i0delta0(2,  6)  =  -1.39016582579643E-04 
i0delta0(l,  6)  =  9.70746619623242E-02 
'  solidity  =1.6 

i0delta0(6,  7)  =  -2.0832321334836  lE-10 
i0delta0(5,  7)  =  3.637575297921 19E-08 
i0delta0(4,  7)  =  -2.51178491095239E-06 
i0delta0(3,  7)  =  7.46391658452694E-05 
i0delta0(2,  7)  = -1.21453189910881E-03 
i0delta0(l,  7)  =  0.11686046293471 
'  solidity  =1.8 

i0delta0(6,  8)  = -1.56143738685847E-10 
i0delta0(5,  8)  =  2.10079058766965E-08 
i0delta0(4,  8)  =  -9.831 16558606056E-07 
i0delta0(3,  8)  =  6.88832778905635E-06 
i0delta0(2,  8)  =  6.30398362773 121E-05 
i0delta0(l,  8)  =  0.124576532031824 
'  solidity  =  2 

i0delta0(6,  9)  =  -3.47201982252689E-10 
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i0delta0(5,  9)  =  6.025 15059199005E-08 
i0delta0(4,  9)  =  -4.22495870910922E-06 
i0delta0(3,  9)  =  1.38861665249124E-04 
i0delta0(2,  9)  =  -2.50961937126704E-03 
i0delta0(l,  9)  =  0.157387592875239 
k  =  9 

ReDim  curve(k) 

ReDim  Xvar(k) 

ReDim  Yvar(k) 
p  =  0.4 

Eor  o  =  1  To  k 

Yvar(o)  =  CurveEit2(i0delta0(  1 ,  o),  i0delta0(2,  o),  i0delta0(3,  o),  i0delta0(4,  o), 
i0delta0(5,  o),  i0delta0(6,  o),  beta(i,  j)) 

Xvar(o)  =  p 

p  =  p  +  0.2 

Next  o 
QuadCoeff  k 

i0d0(i,  j)  =  Interp(sigma(i,  j)) 

If  slopegraph  =  True  Then 

'  NACA  65-(A10)-series  blades  as  equivalent  circular  arc 
'  Load  constants 
'  solidity  =  0.4 

slopemn(7,  1)  =  2.0833333323638  IE- 1 1 
slopemn(6,  1)  =  4.3 108974340537  lE-09 
slopemn(5,  1)  =  3.42948717823449E-07 
slopemn(4,  1)  =  -1.2816870624599E-05 
slopemn(3,  1)  =  1.47926864929815E-04 
slopemn(2,  1)  =  -3.97852565922108E-03 
slopemn(l,  1)  =  0.53497377639842 
'  solidity  =  0.6 

slopemn(7,  2)  =  -3.47222221982991E-12 
slopemn(6,  2)  =  7.772435880668E-10 
slopemn(5,  2)  =  -5.87606836521815E-08 
slopemn(4,  2)  =  1.3414189909966E-06 
slopemn(3,  2)  =  -5.4709838117617E-05 
slopemn(2,  2)  =  -2.55785261049368E-03 
slopemn(l,  2)  =  0.675027681002291 
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'  solidity  =  0.8 

slopemn(7,  3)  =  -5.0821976835258E-21 
slopemn(6,  3)  =  -1.28205 128 154327E-10 
slopemn(5,  3)  =  1.20192307773159E-08 
slopemn(4,  3)  =  -6.92016318915023E-07 
slopemn(3,  3)  =  -3.3 10751 65860284E-05 
slopemn(2,  3)  =  -2.12211544896945E-03 
slopemn(l,  3)  =  0.750014569226494 
'  solidity  =1.0 

slopemn(7,  4)  =  -6.94444444304795E-12 
slopemn(6,  4)  =  1.20192307712444E-09 
slopemn(5,  4)  =  -7.6655982672591E-08 
slopemn(4,  4)  =  1.4852855381 1644E-06 
slopemn(3,  4)  =  4.01670539815768E-05 
slopemn(2,  4)  =  -1.64423084743248E-03 
slopemn(l,  4)  =  0.794953380497141 
'  solidity  =1.2 

slopemn(7,  5)  =  -2.77777777637215E-1 1 
slopemn(6,  5)  =  5.48076922735063E-09 
slopemn(5,  5)  =  4.03579059504722E-07 
slopemn(4,  5)  =  1.270396268893 14E-05 
slopemn(3,  5)  =  -2.01934244898894E-04 
slopemn(2,  5)  =  -3.7852573041 1639E-04 
slopemn(l,  5)  =  0.824973776832628 
'  solidity  =1.4 

slopemn(7,  6)  =  6.94444443457762E-12 
slopemn(6,  6)  =  -9.1346153644617E-10 
slopemn(5,  6)  =  1.5758547106115E-08 
slopemn(4,  6)  =  1.08027389256193E-06 
slopemn(3,  6)  =  -7.76942489579824E-05 
slopemn(2,  6)  =  -3.57051379069162E-04 
slopemn(l,  6)  =  0.850023310682502 
'  solidity  =1.6 

slopemn(7,  7)  =  -1.38888888844018E-1 1 
slopemn(6,  7)  =  3.07692307570384E-09 
slopemn(5,  7)  =  -2.60683760439084E-07 
slopemn(4,  7)  =  9.44930068769168E-06 
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slopemn(3,  7)  =  -1.87603922285007E-04 
slopemn(2,  7)  =  6.39102459899732E-04 
slopemn(l,  7)  =  0.870029138228418 
'  solidity  =1.8 

slopemn(7,  8)  =  -6.94444444643608E-12 
slopemn(6,  8)  =  1.8429487181 1291E-09 
slopemn(5,  8)  =  -1.7841880317615E-07 
slopemn(4,  8)  =  6.76354894046938E-06 
slopemn(3,  8)  =  -1.36371890420151E-04 
slopemn(2,  8)  =  3.83012708653041E-04 
slopemn(l,  8)  =  0.889956294448098 
'  solidity  =  2 

slopemn(7,  9)  =  -1.04166666679601E-1 1 
slopemn(6,  9)  =  2.363782051 1721E-09 
slopemn(5,  9)  =  -2.13141025329211E-07 
slopemn(4,  9)  =  8.38760196941735E-06 
slopemn(3,  9)  =  -1.82296764080547E-04 
slopemn(2,  9)  =  1.1642627027868E-03 
slopemn(l,  9)  =  0.900013112661114 
k  =  9 

ReDim  curve(k) 

ReDim  Xvar(k) 

ReDim  Yvar(k) 
p=  0.4 

Eor  o  =  1  To  k 

Yvar(o)  =  CurveEitl(slopemn(l,  o),  slopemn(2,  o),  slopemn(3,  o), 
slopemn(4,  o),  slopemn(5,  o),  slopemn(6,  o),  slopemn(7,  o), 
beta(i,  j)) 

Xvar(o)  =  p 

p  =  p  +  0.2 

Next  o 
QuadCoeff  k 

onemn(i,  j)  =  Interp(sigma(i,  j)) 

Else 

'  Circular  arc  mean  line  blades 
'  Load  constants 
'  solidity  =  0.4 

slopemn(7,  1)  =  4. 1666666643941  lE-11 

71 


slopemn(6,  1)  =  8.14102563624798E-09 
slopemn(5,  1)  =  -5.809294868 14364E-07 
slopemn(4,  1)  =  1.81949300572665E-05 
slopemn(3,  1)  =  -2.4787 150325345E-04 
slopemn(2,  1)  =  -3.339423083 14938E-03 
slopemn(l,  1)  =  0.350002913822904 
'  solidity  =  0.6 

slopemn(7,  2)  =  6.94444443965982E-12 
slopemn(6,  2)  =  -1.4903846 14550 lE-09 
slopemn(5,  2)  =  1.06303418803688E-07 
slopemn(4,  2)  =  -3.08493589651349E-06 
slopemn(3,  2)  =  -5.89209358992093E-06 
slopemn(2,  2)  =  -2.75032054841517E-03 
slopemn(l,  2)  =  0.545032051554482 
'  solidity  =  0.8 

slopemn(7,  3)  =  -1.04166666603368E-1 1 
slopemn(6,  3)  =  1.56249999787653E-09 
slopemn(5,  3)  =  -7.29166664409364E-08 
slopemn(4,  3)  =  1.18371202617595E-07 
slopemn(3,  3)  =  2. 3873 106883343  lE-05 
slopemn(2,  3)  =  -3.17708338894818E-03 
slopemn(l,  3)  =  0.649981061002052 
'  solidity  =1.0 

slopemn(7,  4)  =  -6.94444443796575E-12 
slopemn(6,  4)  =  7.21153844973621E-10 
slopemn(5,  4)  =  -2.93803398276893E-09 
slopemn(4,  4)  =  -2.18458625589335E-06 
slopemn(3,  4)  =  3.707653560525 16E-05 
slopemn(2,  4)  =  -2.16987186519191E-03 
slopemn(l,  4)  =  0.715017483006385 
'  solidity  =1.2 

slopemn(7,  5)  =  3.47222221644178E-12 
slopemn(6,  5)  =  -1.3221 1538505372E-09 
slopemn(5,  5)  =  1. 4369658 127289E-07 
slopemn(4,  5)  =  -6.95767774239 187E-06 
slopemn(3,  5)  =  1.063561 10350703E-04 
slopemn(2,  5)  =  -2.01842956960263E-03 
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slopemn(l,  5)  =  0.760024767461402 
'  solidity  =1.4 

slopemn(7,  6)  =  -6.7762635780344E-21 
slopemn(6,  6)  =  -1.2820512967221E-10 
slopemn(5,  6)  =  1.20192308605827E-08 
slopemn(4,  6)  =  -8.4353 1474004067E-07 
slopemn(3,  6)  =  -1.71984251977619E-05 
slopemn(2,  6)  =  -7.05448809185327E-04 
slopemn(l,  6)  =  0.794938811805281 
'  solidity  =1.6 

slopemn(7,  7)  =  1.73611111042317E-11 
slopemn(6,  7)  =  -3.34 1346 15 158702E-09 
slopemn(5,  7)  =  2.297008545038E-07 
slopemn(4,  7)  =  -7.7196241337063E-06 
slopemn(3,  7)  =  9.70923186827122E-05 
slopemn(2,  7)  =  -1.42612189358715E-03 
slopemn(l,  7)  =  0.825036422577128 
'  solidity  =1.8 

slopemn(7,  8)  =  1.7361 1110991495E- 11 
slopemn(6,  8)  =  -3.37339743378823E-09 
slopemn(5,  8)  =  2.35309829010877E-07 
slopemn(4,  8)  =  -7.93524184317107E-06 
slopemn(3,  8)  =  9.2248204595080 lE-05 
slopemn(2,  8)  =  -8.26442413767836E-04 
slopemn(l,  8)  =  0.844992716324796 
'  solidity  =  2 

slopemn(7,  9)  =  3.47222221 135958E- 12 
slopemn(6,  9)  =  -5.84935897461614E-10 
slopemn(5,  9)  =  3.55235043048019E-08 
slopemn(4,  9)  =  -2.036349075 14072E-06 
slopemn(3,  9)  =  3.84146771352789E-05 
slopemn(2,  9)  =  -6.8605780506914E-04 
slopemn(l,  9)  =  0.859969406337427 
k  =  9 

ReDim  curve(k) 

ReDim  Xvar(k) 

ReDim  Yvar(k) 
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p  =  0.4 

For  o  =  1  To  k 

Yvar(o)  =  CurveFitl(slopemn(l,  o),  slopemn(2,  o),  slopemn(3,  o), 
slopemn(4,  o),  slopemn(5,  o),  slopemn(6,  o),  slopemn(7,  o), 
beta(i,  j)) 

Xvar(o)  =  p 

p  =  p  +  0.2 

Next  o 
QuadCoeff  k 

onemn(i,  j)  =  Interp(sigma(i,  j)) 

End  If 

'  Variation  of  thickness-correction  factor  Kt  for  camber  calculation 
'  Load  constants 
'  betal  =  10 

Kit(7,  1)  =  816993.464355469 
Kit(6,  1)=  -351150.076196289 
Kit(5,  1)  =  61737.9966625977 
Kit(4,  1)=  -5315.17809592133 
Kit(3,  1)  =  142.41942977427 
Kit(2,  1)=  13.4581061812941 
Kit(L  1)  =  6.500000 14506098E-02 
'  betal  =  20 

Kit(7,  2)  =  1077614.37939453 
Kit(6,  2)  =  -459489.065126953 
Kit(5,  2)  =  79133.1385789185 
Kit(4,  2)  =  -6710.73264049728 
Kit(3,  2)  =  201.572413791658 
Kit(2,  2)  =  12.2108848959033 
Kit(l,  2)  =  7.74781832012336E-02 
'  betal  =  30 

Kit(7,  3)  =  796568.628417969 
Kit(6,  3)  =  -338565.234250488 
Kit(5,  3)  =  58401.6780760498 
Kit(4,  3)=  -4855.18434412186 
Kit(3,  3)  =  103.679839050816 
Kit(2,  3)  =  15.0572606819245 
Kit(l,  3)  =  6.05909105273496E-02 
'  betal  =  40 
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Kit(7,  4)  =  -183823.529541016 
Kit(6,  4)  =  137302.036435547 
Kit(5,4)=  -32831.3536986084 
Kit(4,  4)  =  3879.52703175537 
Kit(3,  4)  =  -335.438802558397 
Kit(2,  4)  =  26.0388210932938 
Kit(l,4)=  -3.53636350253408E-02 
'  betal  =  50 

Kit(7,  5)  =  -449346.408691406 
Kit(6,  5)  =  260840.875805664 
Kit(5,  5)  =  -56572.084102478 
Kit(4,  5)  =  6297.74699841156 
Kit(3,  5)  =  -479.565571368521 
Kit(2,  5)  =  30.7062532322958 
Kit(l,  5)  =  -7.44999987640929E-02 
'  betal  =  60 

Kit(7,  6)  =  2879901.96362305 
Kit(6,  6)=  -1181513.95310303 
Kit(5,  6)  =  193177.319282959 
Kit(4,  6)  =  -15656.7577867749 
Kit(3,  6)  =  525.58354015793 
Kit(2,  6)  =  8.91909924834078 
Kit(l,  6)  =  0.126909092425712 
'  betal  =  70 

Kit(7,  7)  =  5187908.49780273 
Kit(6,  7)  =  -1994626.6972876 
Kit(5,  7)  =  298546.694419434 
Kit(4,  7)  =  -21587.099794791 
Kit(3,  7)  =  600.425619817768 
Kit(2,  7)  =  13.8798799771508 
Kit(l,  7)  =  9.83181832073486E-02 
k  =  7 

ReDim  curve(k) 

ReDim  Xvar(k) 

ReDim  Yvar(k) 

p=  10 

Eor  o  =  1  To  k 
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Yvar(o)  =  CurveFitl(Kit(l,  o),  Kit(2,  o),  Kit(3,  o),  Kit(4,  o),  Kit(5,  o),  Kit(6,  o), 
Kit(7,  o),  tc(i,  j)) 

Xvar(o)  =  p 

p  =  p  +  10 

Next  o 

QuadCoeff  k 

ktbar(i,  j)  =  Interp(beta(i,  j)) 

camber(i,  j)  =  camberfunc(beta(i,  j),  beta((i  -  1),  j),  ksh,  ktbar(i,  j),  i0d0(i,  j),  onemn(i,  j)) 

'  Zero -camber  incidence  angle 
'  Load  constants 
'  solidity  =  0.4 

i0(6,  1)  =  -4.04805786872846E-14 
i0(5,  1)  =  3.30321577382553E-10 
i0(4,  1)  =  -1.6 119368995 1489E-07 
i0(3,  1)=  1.602666840 128  llE-05 
i0(2,  1)  =  -6.12705630658184E-04 
i0(l,  1)  =  3.964147943543 14E-02 
'  solidity  =  0.6 

i0(6,  2)  =  2.77737297393952E-10 
i0(5,  2)  =  -5.800301 1805044E-08 
i0(4,  2)  =  4.49158409132622E-06 
i0(3,  2)  =  -1.58973331821244E-04 
i0(2,  2)  =  2.41923882140327E-03 
i0(l,  2)  =  3.772481264741 15E-02 
'  solidity  =  0.8 

i0(6,  3)  =  -2.77737297421058E-10 
i0(5,  3)  =  5.38363451353663E-08 
i0(4,  3)  =  -3.86658408690899E-06 
i0(3,  3)  =  1.2355666447661E-04 
i0(2,  3)  =  -1.73173878386024E-03 
i0(l,  3)  =  7.33585 199368463E-02 
'  solidity  =1.0 

i0(6,  4)  =  -6.24048706766443E-10 
i0(5,  4)  =  1. 33904 109689276E-07 
i0(4,  4)  =  -1.07952815886492E-05 
i0(3,  4)  =  3.96289954153417E-04 
i0(2,  4)  =  -6.38683407 169083E-03 
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iO(l,  4)  =  0.11471689445807 
'  solidity  =1.2 

i0(6,  5)  =  1.3907105163107E-10 
i0(5,  5)  =  -2.856978047 1943 lE-08 
i0(4,  5)  =  2.01009383982154E-06 
i0(3,  5)  =  -5.64950093888683E-05 
i0(2,  5)  =  3.77314309162102E-04 
i0(l,  5)  =  0.100321674263796 
'  solidity  =1.4 

i0(6,  6)  =  -1.0411606600198  IE- 10 
i0(5,  6)  =  2.042043 13875023E-08 
i0(4,  6)  =  -1.67350788249365E-06 
i0(3,  6)  =  6.850833028693 16E-05 
i0(2,  6)  =  - 1.4205762442998  lE-03 
i0(l,  6)  =  0.123885649683871 
'  solidity  =1.6 

i0(6,  7)  =  -6.59874024969446E-10 
i0(5,  7)  =  1.26238706023296E-07 
i0(4,  7)  =  -9.07669856742288E-06 
i0(3,  7)  =  3.00725006795233E-04 
i0(2,  7)  =  -4.56882667 162972E-03 
i0(l,  7)  =  0.153343047833914 
'  solidity  =1.8 

i0(6,  8)  =  -2.42843032737926E-10 
i0(5,  8)  =  3.9932478421 1563E-08 
i0(4,  8)  =  -2.4176220057015E-06 
i0(3,  8)  =  6.48183223432852E-05 
i0(2,  8)  =  -8.29475889986497E-04 
i0(l,  8)  =  0.149069730896372 
'  solidity  =  2 

i0(6,  9)  =  -6.5928705639891E-10 
i0(5,  9)  =  1.27699043138418E-07 
i0(4,  9)  =  -9.55189003892798E-06 
i0(3,  9)  =  3.42296645044371E-04 
i0(2,  9)  =  -5.90334486059874E-03 
i0(l,  9)  =  0.198583259410952 
k  =  9 
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ReDim  curve(k) 

ReDim  Xvar(k) 

ReDim  Yvar(k) 
p  =  0.4 

For  o  =  1  To  k 

Yvar(o)  =  CurveFit2(iO(l,  o),  i0(2,  o),  i0(3,  o),  i0(4,  o),  i0(5,  o),  i0(6,  o),  beta(i, 

j)) 

Xvar(o)  =  p 

p  =  p  +  0.2 

Nexto 
QuadCoeff  k 

i010(i,  j)  =  Interp(sigma(i,  j)) 

'  Minimum  loss  incidence  angle  slope  factor 
'  Load  constants 
'  solidity  =  0.4 

slopen(7,  1)  =  5.20833333567409E-12 
slopen(6,  1)  =  -1.00560897412991E-09 
slopen(5,  1)  =  7.05128204714356E-08 
slopen(4,  1)  =  -2.35449445540326E-06 
slopen(3,  1)  =  7.51638881979488E-06 
slopen(2,  1)  = -3.64078518850874E-03 
slopen(l,  1)  =  -5.00123838208282E-02 
'  solidity  =  0.6 

slopen(7,  2)  =  -1.7361 11 10737386E- 12 
slopen(6,  2)  =  2.604 166663 12754E- 10 
slopen(5,  2)  =  -2.2569444452225  lE-08 
slopen(4,  2)  =  1.0866477326843E-06 
slopen(3,  2)  =  -5.59643317217251E-05 
slopen(2,  2)  =  -2.28645827399987E-03 
slopen(l,  2)  =  -4.49905306238207E-02 
'  solidity  =  0.8 

slopen(7,  3)  =  6.94444444474202E-12 
slopen(6,  3)  =  -1.85897435877984E-09 
slopen(5,  3)  =  1.70806623789321E-07 
slopen(4,  3)  = -7.15544871 177087E-06 
slopen(3,  3)  =  1.05446046163138E-04 
slopen(2,  3)  =  -2.7290063578 107E-03 
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slopen(l,  3)  =  -4.00160259214886E-02 
'  solidity  =1.0 

slopen(7,  4)  =  -5.2083333280508E-12 
slopen(6,  4)  =  9.57532051044929E-10 
slopen(5, 4)  =  -7.77243589489274E-08 
slopen(4,  4)  =  3.08220425893069E-06 
slopen(3,  4)  =  -8.76518800509984E-05 
slopen(2,  4)  =  -8.65945466387075E-04 
slopen(l,  4)  =  -3. 4996358054 1726E-02 
'  solidity  =1.2 

slopen(7,  5)  =  2.60416666588887E-11 
slopen(6,  5)  =  -5.733 1730743057 lE-09 
slopen(5,  5)  =  4.603365382355 17E-07 
slopen(4,  5)  =  -1.68285620523 179E-05 
slopen(3,  5)  =  2.388239357 17392E-04 
slopen(2,  5)  =  -2.12556085602955E-03 
slopen(l,  5)  =  -2.99817892565741E-02 
'  solidity  =1.4 

slopen(7,  6)  =  1.56249999943168E-11 
slopen(6,  6)  =  -3.43349358927075E-09 
slopen(5,  6)  =  2.74038461367532E-07 
slopen(4,  6)  =  -1.01127258078648E-05 
slopen(3,  6)  =  1.33344623499454E-04 
slopen(2,  6)  =  -1.06402239998715E-03 
slopen(l,  6)  =  -2.49992717449814E-02 
'  solidity  =1.6 

slopen(7,  7)  =  1.56249999951638E-11 
slopen(6,  7)  =  -3.32131410091599E-09 
slopen(5,  7)  =  2.54407051 133998E-07 
slopen(4,  7)  =  -8.82776077215652E-06 
slopen(3,  7)  =  9.46172051285998E-05 
slopen(2,  7)  =  -2.87900609066583E-04 
slopen(l,  7)  =  -1.99978148450413E-02 
'  solidity  =1.8 

slopen(7,  8)  =  -5.20833332889783E-12 
slopen(6,  8)  =  1.11778846096679E-09 
slopen(5,  8)  =  -9.01442307604805E-08 
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slopen(4,  8)  =  2.80612616787579E-06 
slopen(3,  8)  =  -5.87438087080727E-05 
slopen(2,  8)  =  3.75240412154199E-04 
slopen(l,  8)  =  -0.01498615981 1361 
'  solidity  =  2 

slopen(7,  9)  =  1.21527777732163E-11 
slopen(6,  9)  = -2.31971 153754913E-09 
slopen(5,  9)  =  1.6159188022391E-07 
slopen(4,  9)  =  -5.68345716178698E-06 
slopen(3,  9)  =  7.840896905975 llE-05 
slopen(2,  9)  =  -3.66426257727426E-04 
slopen(l,  9)  =  -9.98615979092321E-03 
k  =  9 

ReDim  curve(k) 

ReDim  Xvar(k) 

ReDim  Yvar(k) 
p  =  0.4 

Eor  o  =  1  To  k 

Yvar(o)  =  CurveEitl(slopen(l,  o),  slopen(2,  o),  slopen(3,  o),  slopen(4,  o), 
slopen(5,  o),  slopen(6,  o),  slopen(7,  o),  beta(i,  j)) 

Xvar(o)  =  p 

p  =  p  -I-  0.2 

Next  o 

QuadCoeff  k 

N(i,  j)  =  Interp(sigma(i,  j)) 

'  Maximum  thickness  correction  factor 
'  Load  constants 
ikit(6)  =-748795.365783691 
ikit(5)  =  243951.67288208 
ikit(4)  =  -28087.2979736328 
ikit(3)=  1612.26135325431 
ikit(2)  =-137.429116554558 
ikit(l)=  18.8187430176185 

ikt(i,  j)  =  CurveEit2(ikit(l),  ikit(2),  ikit(3),  ikit(4),  ikit(5),  ikit(6),  tc(i,  j)) 
i0ref(i,  j)  =  ksh  *  ikt(i,  j)  *  i010(i,  j) 
i2d  =  i0ref(i,  j)  -i-  N(i,  j)  *  camber(i,  j) 

If  i  =  2  Then 
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If  ksh  =  0.7  Then 

icor(i,  j)  =  0.7238  *  Mw(i  -  1,  j)  +  7.5481  +  i2d 
Elself  ksh  =1.1  Then 

icor(i,j)=  1.3026  *Mw(i-  1,  j)  +  5.738  +  i2d 
Elself  ksh  =  1  Then 

icor(i,  j)  =  i2d 

End  If 

Elself  i  =  3  Then 

If  ksh  =  0.7  Then 

icor(i,  j)  =  0.7238  *  M(i  -  1,  j)  +  7.5481  +  i2d 
Elself  ksh  =1.1  Then 

icor(i,  j)  =  1.3026  *  M(i  -  1,  j)  +  5.738  +  i2d 
Elself  ksh  =  1  Then 

icor(i,  j)  =  i2d 

End  If 

End  If 

'  Zero -camber  deviation  angle 
'  Load  constants 
'  solidity  =  0.4 

d0(6,  1)  =  -6.95861265439764E-11 
d0(5,  1)  =  1.57394588647108E-08 
d0(4,  1)=  -1.36278903151155E-06 
d0(3,  1)  =  5.7135006738207E-05 
d0(2,  1)  =  - 1.0639 14 18466956E-03 
d0(l,  1)  =  1.59951786608872E-02 
'  solidity  =  0.6 

d0(6,  2)  =  -7.08410313374729E-14 
d0(5,  2)  =  -4.63603904397869E-10 
d0(4,  2)  =  1.86661040402214E-07 
d0(3,  2)  =  -1.62241633923088E-05 
d0(2,  2)  =  6.30890134971196E-04 
d0(l,  2)  =  5.06008922957335E-03 
'  solidity  =  0.8 

d0(6,  3)  =  -1.72852100539214E-11 
d0(5,  3)  =  3.54731371 186856E-09 
d0(4,  3)  =  -1.83872938541718E-07 
d0(3,  3)  =  1.7208295339799E-06 
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d0(2,  3)  =  2.60107809708643E-04 
dO(l,  3)  =  9.9534754635897 lE-03 
'  solidity  =1.0 

d0(6,  4)  =  1.91 1 898054967  04E- 10 
d0(5,  4)  =  -3.92754785100841E-08 
d0(4,  4)  =  3.08863831766093E-06 
d0(3,4)=  -1. 106225 10739233E-04 
d0(2,  4)  =  2.00666089165225E-03 
d0(l,  4)  =  9.999633 12096952E-04 
'  solidity  =1.2 

d0(6,  5)  =  -2.95336231991993E-10 
d0(5,  5)  =  6.09853 177350322E-08 
d0(4, 5)  =  -4.62262480382947E-06 
d0(3,  5)  =  1.62505009349445E-04 
d0(2,  5)  =  -2.2878080447981  lE-03 
d0(l,  5)  =  2.54709629980425E-02 
'  solidity  =1.4 

d0(6,  6)  =  -1.75483338526582E-11 
d0(5,  6)  =  1.5277372943101  lE-09 
d0(4,  6)  =  2.26701400096729E-07 
d0(3,  6)  =  -2.11891579624535E-05 
d0(2,  6)  =  8.81687814285215E-04 
d0(l,  6)  =  7.87309303404982E-03 
'  solidity  =1.6 

d0(6,  7)  =  -3.46716215172916E-11 
d0(5,  7)  =  7.92043135938725E-09 
d0(4,  7)  =  -5.62396773595708E-07 
d0(3,  7)  =  2.2674997246952E-05 
d0(2,  7)  =  -1.57381818780777E-04 
d0(l,  7)  =  0.017343983408864 
'  solidity  =1.8 

d0(6,  8)  =  6.30485119950015E-11 
d0(5,  8)  =  -1.41425240546278E-08 
d0(4,  8)  =  1.37167452307629E-06 
d0(3,  8)  =  -5.62863600421792E-05 
d0(2,  8)  =  1.32778646161569E-03 
d0(l,  8)  =  8.467 11704980407E-03 


82 


'  solidity  =  2 

d0(6,  9)  =  8.67296545047297E-1 1 
d0(5,  9)  =  -1.60473137209016E-08 
d0(4,  9)  =  1.29498405332384E-06 
d0(3,  9)  =  -4.338749676691 12E-05 
d0(2,  9)  =  8.78086666489253E-04 
dO(l,  9)  =  1.43381905636488E-02 
k  =  9 

ReDim  curve(k) 

ReDim  Xvar(k) 

ReDim  Yvar(k) 
p  =  0.4 

Eoro  =  1  To  k 

Yvar(o)  =  CurveEit2(dO(l,  o),  d0(2,  o),  d0(3,  o),  d0(4,  o),  d0(5,  o),  d0(6,  o), 
beta(i,  j)) 

Xvar(o)  =  p 

p  =  p  -I-  0.2 

Next  o 
QuadCoeff  k 

d010(i,  j)  =  Interp(sigma(i,  j)) 

If  slopegraph  =  True  Then 

'  NACA  65-(A10)-series  blades  as  equivalent  circular  arc 
'  Load  constants 
'  solidity  =  0.4 

slopem(7,  1)  =  6.94444443457762E-12 
slopem(6,  1)  = -1.49038461433326E-09 
slopem(5,  1)  =  1.06303418845322E-07 
slopem(4,  1)  = -3.16069348027526E-06 
slopem(3,  1)  =  7.7062452689347  lE-05 
slopem(2,  1)  = -7.91987265529315E-04 
slopem(l,  1)  =  0.412494173012831 
'  solidity  =  0.6 

slopem(7,  2)  =  2.08333333206735E-1 1 
slopem(6,  2)  = -4.48717948388355E-09 
slopem(5,  2)  =  3.5296474354407  lE-07 
slopem(4,  2)  = -1.22224650365 155E-05 
slopem(3,  2)  =  2.05185752662373E-04 
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slopem(2,  2)  =  -8.88621860667627E-04 
slopem(l,  2)  =  0.277498543514405 
'  solidity  =  0.8 

slopem(7,  3)  =  1.0416666661 1838E-1 1 
slopem(6,  3)  = -2.20352563908 184E-09 
slopem(5,  3)  =  1.64262820512295E-07 
slopem(4,  3)  = -4.7716346145421  lE-06 
slopem(3,  3)  =  6.08733981977139E-05 
slopem(2,  3)  =  2.12339688573593E-04 
slopem(l,  3)  =  0.210016025961522 
'  solidity  =1.0 

slopem(7,  4)  =  3.47222221813584E-12 
slopem(6,  4)  = -1.19391025581519E-09 
slopem(5,  4)  =  1.31677350481696E-07 
slopem(4,  4)  =  -5.84899476141487E-06 
slopem(3,  4)  =  1.26963627380405E-04 
slopem(2,  4)  =  -5.62980817107928E-04 
slopem(l,  4)  =  0.170010198410182 
'  solidity  =1.2 

slopem(7,  5)  =  5.2083333305919E-12 
slopem(6,  5)  = -1.5 184294867472 lE-09 
slopem(5,  5)  =  1. 550480769885 18E-07 
slopem(4,  5)  =  -6.59036276573488E-06 
slopem(3,  5)  =  1.3337230542021E-04 
slopem(2,  5)  =  -5.56330170923047E-04 
slopem(l,  5)  =  0.142489073668344 
'  solidity  =1.4 

slopem(7,  6)  =  2.256944443 14355E-1 1 
slopem(6,  6)  =  -4.79567307436549E-09 
slopem(5,  6)  =  3.83947649437721E-07 
slopem(4,  6)  =  -1.38976908479194E-05 
slopem(3,  6)  =  2.37974820819886E-04 
slopem(2,  6)  = -1.0776442692304E-03 
slopem(l,  6)  =  0.122499271780747 
'  solidity  =1.6 

slopem(7,  7)  = -1.04166666620308E-11 
slopem(6,  7)  =  1.94711538320687E-09 
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slopem(5,  7)  = -1.29807692192396E-07 
slopem(4,  7)  =  4.01260197557818E-06 
slopem(3,  7)  =  -4.06300984678865E-05 
slopem(2,  7)  =  3.55929451 188786E-04 
slopem(l,  7)  =  0.107513112087279 
'  solidity  =1.8 

slopem(7,  8)  =  1.21527777736398E-1 1 
slopem(6,  8)  =  -2.80048076742312E-09 
slopem(5,  8)  =  2.45726495706755E-07 
slopem(4,  8)  =  -9.71317744902223E-06 
slopem(3,  8)  =  1.8301998394321  lE-04 
slopem(2,  8)  = -8.29567340645099E-04 
slopem(l,  8)  =  9.50021854990979E-02 
'  solidity  =  2 

slopem(7,  9)  = -3.47222222152398E-12 
slopem(6,  9)  =  6.009615380201 19E-10 
slopem(5,  9)  = -3.31 196580577453E-08 
slopem(4,  9)  =  8.84688224900287E-07 
slopem(3,  9)  =  4.18366859378239E-07 
slopem(2,  9)  =  1.88301250858558E-04 
slopem(l,  9)  =  8.74956295422606E-02 
k  =  9 

ReDim  curve(k) 

ReDim  Xvar(k) 

ReDim  Yvar(k) 
p  =  0.4 

Eor  o  =  1  To  k 

Yvar(o)  =  CurveEitl(slopem(l,  o),  slopem(2,  o),  slopem(3,  o), 
slopem(4,  o),  slopem(5,  o),  slopem(6,  o),  slopem(7,  o),  beta(i, 

j)) 

Xvar(o)  =  p 
p  =  p  -I-  0.2 

Next  o 

QuadCoeff  k 

dm(i,  j)  =  Interp(sigma(i,  j)) 

Else 

'  Circular  arc  mean  line  blades 
'  Load  constants 

85 


'  solidity  =  0.4 

slopem(7,  1)  = -5.0821976835258E-21 
slopem(6,  1)  =  6.73076922756005E-10 
slopem(5,  1)  = -1.1258012813653E-07 
slopem(4,  1)  =  6.10941141054866E-06 
slopem(3,  1)  = -1.2890806971 1997E-04 
slopem(2,  1)  =  1.04214733619301E-03 
slopem(l,  1)  =  0.602489802490481 
'  solidity  =  0.6 

slopem(7,  2)  =  5.20833332212157E-12 
slopem(6,  2)  = -1.0857371781 1505E-09 
slopem(5,  2)  =  8.45352564893265E-08 
slopem(4,  2)  =  -3.08293269490889E-06 
slopem(3,  2)  =  5.90424689903557E-05 
slopem(2,  2)  =  -9.991994829761  lE-05 
slopem(l,  2)  =  0.409991987649079 
'  solidity  =  0.8 

slopem(7,  3)  =  1.73611110229166E-12 
slopem(6,  3)  = -1.9631410169349E-10 
slopem(5,  3)  =  9.34829047505303E-10 
slopem(4,  3)  =  6.51405884788403E-07 
slopem(3,  3)  = -1.73966821961358E-05 
slopem(2,  3)  =  6.07932628 156505E-04 
slopem(l,  3)  =  0.310002185699076 
'  solidity  =1.0 

slopem(7,  4)  =  -6.94444444474202E-12 
slopem(6,  4)  =  1.66666666600625E-09 
slopem(5,  4)  =  -1.47569444275408E-07 
slopem(4,  4)  =  6.10795453948043E-06 
slopem(3,  4)  =  -1.06196337540609E-04 
slopem(2,  4)  =  1.09791661 13998E-03 
slopem(l,  4)  =  0.249981060934445 
'  solidity  =1.2 

slopem(7,  5)  = -1.215277777 19458E- 11 
slopem(6,  5)  =  2.57612179309884E-09 
slopem(5,  5)  = -2.01255341752971E-07 
slopem(4,  5)  =  7.285292825898 18E-06 
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slopem(3,  5)  =  -1.1006325 1923975E-04 
slopem(2,  5)  =  1.03774033499349E-03 
slopem(l,  5)  =  0.210013840616803 
'  solidity  =1.4 

slopem(7,  6)  = -6.94444444135388E-12 
slopem(6,  6)  =  1.53846153871928E-09 
slopem(5,  6)  = -1.199252 13509786E-07 
slopem(4,  6)  =  4.09965034364745E-06 
slopem(3,  6)  = -4.48436278688291E-05 
slopem(2,  6)  =  3.82051237181713E-04 
slopem(l,  6)  =  0.182514569024661 
'  solidity  =1.6 

slopem(7,  7)  = -2.25694444348236E-11 
slopem(6,  7)  =  4.50721 153650614E-09 
slopem(5,  7)  =  -3.28258546761218E-07 
slopem(4,  7)  =  1.07734192780207E-05 
slopem(3,  7)  = -1.40615408163569E-04 
slopem(2,  7)  =  8.60176240280452E-04 
slopem(l,  7)  =  0.160005099306115 
'  solidity  =1.8 

slopem(7,  8)  =  1.7361 1110991495E-12 
slopem(6,  8)  =  -5.48878204605779E-10 
slopem(5,  8)  =  6.78418804189729E-08 
slopem(4,  8)  =  -3.58591929172647E-06 
slopem(3,  8)  =  9.18654094732574E-05 
slopem(2,  8)  =  -3.68509653696947E-04 
slopem(l,  8)  =  0.142494901151593 
'  solidity  =  2 

slopem(7,  9)  =  -5.20833333567409E-12 
slopem(6,  9)  =  1.14983974338484E-09 
slopem(5,  9)  =  -9.054487 17253288E-08 
slopem(4,  9)  =  3.42894084548462E-06 
slopem(3,  9)  = -5.49887086620515E-05 
slopem(2,  9)  =  7.85977528323656E-04 
slopem(l,  9)  =  0.1275007286422 
k  =  9 

ReDim  curve(k) 
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ReDim  Xvar(k) 

ReDim  Yvar(k) 
p  =  0.4 

For  o  =  1  To  k 

Yvar(o)  =  CurveFitl(slopem(l,  o),  slopem(2,  o),  slopem(3,  o), 
slopem(4,  o),  slopem(5,  o),  slopem(6,  o),  slopem(7,  o),  beta(i, 

j)) 

Xvar(o)  =  p 

p  =  p  +  0.2 

Next  o 
QuadCoeff  k 

dm(i,  j)  =  Interp(sigma(i,  j)) 

End  If 

'  Maximum  thickness  correction  factor 
'  Load  constants 
dkit(6)  =  618823.625244141 
dkit(5)  =  -202775.302703857 
dkit(4)  =  25013.8597869873 
dkit(3)=  -1269.01561832427 
dkit(2)  =  41.3428950682282 
dkit(l)  =  7.56794627627824 

dkt(i,  j)  =  CurveFit2(ikit(l),  ikit(2),  ikit(3),  ikit(4),  ikit(5),  ikit(6),  tc(i,  j)) 
d0ref(i,  j)  =  ksh  *  dkt(i,  j)  *  d010(i,  j) 
dref(i,  j)  =  d0ref(i,  j)  -i-  dm(i,  j)  *  camber(i,  j) 

End  Sub 

14.  Module  mdlinterpolation  Code 

Option  Explicit 
Option  Base  1 
Public  AinterpO  As  Double 
Public  BinterpO  As  Double 
Public  CinterpO  As  Double 
Dim  o  As  Integer 

Public  Sub  QuadCoeff(k  As  Integer) 

ReDim  Ainterp(k) 

ReDim  Binterp(k) 

ReDim  Cinterp(k) 

Eor  o  =  2  To  k  -  1 
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Ainterp(o)  =  interpA(Xvar(o  +  1),  Xvar(o  -  1),  Xvar(o),  Yvar(o),  Yvar(o  +  1), 
Yvar(o  -  1)) 

Binterp(o)  =  interpB(Yvar(o),  Yvar(o  -  1),  Xvar(o),  Xvar(o-  1),  Ainterp(o)) 
Cinterp(o)  =  interpC(Yvar(o),  Ainterp(o),  Binterp(o),  Xvar(o)) 

Next  o 

Ainterp(l)  =  Ainterp(2) 

Binterp(l)  =  Binterp(2) 

Cinterp(l)  =  Cinterp(2) 

Ainterp(k)  =  Ainterp(k-  1) 

Binterp(k)  =  Binterp(k  -  1) 

Cinterp(k)  =  Cinterp(k  -  1) 

For  o  =  1  To  k  -  1 

Ainterp(o)  =  (Ainterp(o)  +  Ainterp(o  +  1))  /  2 
Binterp(o)  =  (Binterp(o)  +  Binterp(o  +  1))  /  2 
Cinterp(o)  =  (Cinterp(o)  +  Cinterp(o  +  1))  /  2 

Next  o 

End  Sub 

Public  Function  Interp(Xval  As  Double) 

0  =  2 
Do 

If  Xval  <  Xvar(o)  Then 

Interp  =  Ainterp(o)  *  Xval  ^2  +  Binterp(o)  *  Xval  +  Cinterp(o) 

Exit  Do 

Else 

If  Xval  >  Xvar(o)  And  Xval  <  Xvar(o  +1)  Then 

Interp  =  Ainterp(o)  *  Xval  ''  2  +  Binterp(o)  *  Xval  + 
Cinterp(o) 

Exit  Do 

Else 

0  =  0+1 

End  If 

End  If 

If  o  =  k  -  I  Then 

Interp  =  Ainterp(o)  *  Xval  *  2  +  Binterp(o)  *  Xval  +  Cinterp(o) 

Exit  Do 

End  If 

Loop  Until  o  =  k 


End  Eunction 
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15.  Module  mdlFunctions  Code 


Option  Exp  licit 
Option  Base  1 

Public  Function  Diffusion(phil  As  Double,  betal  As  Double,  beta2  As  Double,  pitch_ml  As 
Double,  solidity  As  Double)  As  Double 

Diffusion  =  1  -  phil  *  (Cos(DegToRad(betal))  /  Cos(DegToRad(beta2)))  + 

(((Tan(DegToRad(betal))  -  pitch_ml  *  phil  *  Tan(DegToRad(beta2)))  * 
Cos(DegToRad(betal)))  /  ((1  +  pitch_ml)  *  solidity)) 

End  Function 

Public  Function  ggl(g  As  Double)  As  Double 
ggl  =  g/(g  -  1) 

End  Function 

Public  Function  gl(g  As  Double)  As  Double 
gl  =  l/(g-  1) 

End  Function 

Public  Function  glg(g  As  Double)  As  Double 
gig  =  (g-  l)/g 
End  Function 

Public  Function  Mach(g  As  Double,  XI  As  Double)  As  Double 
Mach  =  Sqr(((2  /  (g  -  1))  *  (XI  2))  /  (1  -  (XI  2))) 

End  Function 

Public  Function  DegToRad(angle  As  Double)  As  Double 
DegToRad  =  angle  *  ((22  /  7)  /  180) 

End  Function 

Public  Function  RadToDeg(radians  As  Double)  As  Double 
RadToDeg  =  radians  *  (180  /  (22  /  7)) 

End  Function 

Public  Function  Arccos(ratio  As  Double)  As  Double 

Arccos  =  RadToDeg(Atn(-ratio  /  Sqr(-ratio  *  ratio  +  1))  +  2  *  Atn(l)) 

End  Function 

Public  Function  Arcsin(ratio  As  Double)  As  Double 
Arcsin  =  RadToDeg(Atn(ratio  /  Sqr( -ratio  *  ratio  +  1))) 

End  Function 

Public  Function  DofReaction(Xthetal  As  Double,  Xul  As  Double,  Xtheta2  As  Double,  XU2  As 
Double)  As  Double 

DofReaction  =  (1  -  0.5  *  ((Xthetal  /  Xul)  +  (Xtheta2  /  XU2))) 

End  Function 
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Public  Function  SFTC(betal  As  Double,  beta2  As  Double,  solidity  As  Double,  height  As  Double, 
tipgap  As  Double,  spacing  As  Double) 

Dim  CD!  As  Double 

Dim  CL  As  Double 

Dim  betainf  As  Double 

betainf  =  Atn((Tan(DegToRad(betal))  +  Tan(DegToRad(beta2)))  /  2) 

CL  =  (2  /  solidity)  *  (Tan(DegToRad(betal))  -  Tan(DegToRad(beta2)))  * 

Cos(DegToRad(betainf)) 

CD!  =  (0.25  *  (CL  *  2)  *  solidity  *  (tipgap  /  height)  *  (1  /  Cos(DegToRad(beta2))))  + 
(0.04  *  (CL  *  2)  *  solidity  *  (spacing  /  height)) 

SFTC  =  CDi  *  ((Cos(DegToRad(betal))  2)  /  (Cos(DegToRad(betainf))  3))  *  solidity 

End  Function 

Public  Function  taufunc(XU2  As  Double,  Xtheta2  As  Double,  Xul  As  Double,  Xthetal  As 
Double)  As  Double 

taufunc  =  1  +  2  *  ((XU2  *  Xtheta2)  -  (Xul  *  Xthetal)) 

End  Eunction 

Public  Eunction  Xfunc(Xzl  As  Double,  alphal  As  Double)  As  Double 
Xfunc  =  Xzl  /  Cos(DegToRad(alphal)) 

End  Eunction 

Public  Eunction  Xthetafunc(rl  As  Double,  r2  As  Double,  Xthetal  As  Double)  As  Double 
Xthetafunc  =  (rl  /  r2)  *  Xthetal 
End  Eunction 

Public  Eunction  Xufunc(rl  As  Double,  r2  As  Double,  Xul  As  Double)  As  Double 
Xufunc  =  (r2  /  rl)  *  Xul 
End  Eunction 

Public  Eunction  Xwfunc(Xzl  As  Double,  betal  As  Double)  As  Double 
Xwfunc  =  Xzl  /  Cos(DegToRad(betal)) 

End  Eunction 

Public  Eunction  Yfunc(Xl  As  Double,  taul  As  Double)  As  Double 
Yfunc  =  XI  /  Sqr(taul) 

End  Eunction 

Public  Eunction  Machz(Ml  As  Double,  alphal  As  Double)  As  Double 
Machz  =  Ml  *  Cos(DegToRad(alphal)) 

End  Eunction 

Public  Eunction  alphafunc(Xthetal  As  Double,  Xzl  As  Double)  As  Double 
alphafunc  =  RadToDeg(Atn(Xthetal  /  Xzl)) 

End  Eunction 

Public  Eunction  betafunc(Xul  As  Double,  Xthetal  As  Double,  Xzl  As  Double)  As  Double 
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betafunc  =  RadToDeg(Atn((Xul  -  Xthetal)  /  Xzl)) 

End  Function 

Public  Function  rhfuncl(rt  As  Double,  A1  As  Double)  As  Double 
rhfuncl  =  Sqr(rt  2  -  (A1  /  (22  /  7))) 

End  Function 

Public  Function  rhfunc2(rni  As  Double,  rt  As  Double)  As  Double 
rhfunc2  =  2  *  rm  -  rt 
End  Function 

Public  Function  rtfunc(Al  As  Double,  rm  As  Double)  As  Double 
rtfunc  =  (A1  /  (4  *  rm  *  (22  /  7)))  +  rm 
End  Function 

Public  Function  rmfunc(rt  As  Double,  rh  As  Double)  As  Double 
rmfunc  =  (rt  +  rh)  /  2 
End  Function 

Public  Function  rhtfunc(rt  As  Double,  rh  As  Double)  As  Double 
rhtfunc  =  rh  /  rt 
End  Function 

Public  Function  rhtfunc2(h  As  Double,  rm  As  Double)  As  Double 
Dim  hrm  As  Double 
hrm  =  h  /  (2  *  rm) 
rhtfunc2  =  (1  -  hrm)  /  (1  +  hrm) 

End  Function 

Public  Function  CurveFitl(constO  As  Double,  const  1  As  Double,  const2  As  Double,  const3  As 
Double,  const4  As  Double,  constS  As  Double,  constb  As  Double,  polyvar  As  Double) 

'  No  y  intercept 

CurveFitl  =  (constb  *  (polyvar  6))  +  (constS  *  (polyvar  5))  +  (const4  *  (polyvar  4)) 
+  (constS  *  (polyvar  3))  +  (const2  *  (polyvar  *  2))  +  (const  1  *  polyvar)  + 
constO 

End  Function 

Public  Function  CurveFit2(constl  As  Double,  const2  As  Double,  constS  As  Double,  const4  As 
Double,  constS  As  Double,  constb  As  Double,  polyvar  As  Double) 

'  y  intercept 

CurveFit2  =  (constb  *  (polyvar  *  6))  +  (constS  *  (polyvar  *  S))  +  (const4  *  (polyvar  4)) 
+  (constS  *  (polyvar  3))  +  (const2  *  (polyvar  *  2))  +  (const  1  *  polyvar) 

End  Function 

Public  Function  camberfunc(beta2  As  Double,  betal  As  Double,  kshape  As  Single,  kthick  As 
Double,  idelta  As  Double,  slope  As  Double) 

camberfunc  =  ((beta2  -  betal)  -  (kshape  *  kthick  *  idelta))  /  slope 

End  Function 
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Public  Function  interpA(Xipl  As  Double,  Ximl  As  Double,  Xi  As  Double,  Yi  As  Double,  Yipl 
As  Double,  Yiml  As  Double) 

interpA  =  (1  /  (Xipl  -  Ximl))  *  (((Yipl  -  Yi)  /  (Xipl  -  Xi))  -  ((Yi  -  Yiml)  /  (Xi  -  Ximl))) 
End  Function 

Public  Function  interpB(Yi  As  Double,  Yiml  As  Double,  Xi  As  Double,  Ximl  As  Double,  Ai  As 
Double) 

interpB  =  ((Yi  -  Yiml)  /  (Xi  -  Ximl))  -  Ai  *  (Xi  +  Ximl) 

End  Function 

Public  Function  interpC(Yi  As  Double,  Ai  As  Double,  Bi  As  Double,  Xi  As  Double) 
interpC  =  Yi  -  Ai  *  Xi  2  -  Bi  *  Xi 
End  Function 

Public  Function  ShockLoss(M  As  Double,  g  As  Double) 

Dim  ystar  As  Double 

ystar  =  (1  /  (4  *  g  *  M  2))  *  (((g  +  1)  *  M  2)  -  (3  -  g)  +  Sqr((g  +  1)  *  (((g  +  1)  *  M  4) 
-  (2  *  (3  -  g)  *  M  2)  +  g  +  9))) 

ShockLoss  =  ((g  +  1)  /  (2  *  g  *  M  *  2  *  ystar  -  (g  -  1)))  gl(g)  *  (((g  +  1)  *  M  2  * 
ystar)  /  (2  +  (g- 1)  *  M  *  2  *  ystar))  *  ggl(g) 

End  Function 

Public  Function  DiffB  (sigma  As  Double,  R  As  Double) 

DiffB  =  2  *  sigma  *  R  *  ( 1  +  R) 

End  Function 

Public  Function  DiffA(sigma  As  Double,  R  As  Double,  D  As  Double,  phi  As  Double,  betal  As 
Double) 

DiffA  =  (1  +  R)  *  sigma  *  ((1  -  D)  /  (phi  *  Cos(DegToRad(betal))))  + 
(Tan(DegToRad(betal))  /  phi) 

End  Function 

Public  Function  SinB2(A  As  Double,  B  As  Double,  R  As  Double) 

SinB2  =  (-B  +  Sqr(B  ''2  +  (4*(A''2-(B''2/(4*R^  2)))  *  (R  2  +  A  2))))  /  (2  *  (R 
2  +  A  2)) 

End  Function 
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APPENDIX  C.  CALCULATION  RESULTS 

Cl.  HAND  CALCULATIONS 


Inputs  derived  from  Sanger  [Ref.  13]. 
rh  =  17.091bm/sec 
Pj  =14.69psia 
T,  =  517.8"R 

Mwu=l-28 
(0=  2836.33rad/sec 
R  =  53.35ft  Ibsf/lbm  -R 
7  =  1.4 


a,.  =  0- 


U|  =  CO  =  2836.33  •  (5.5  inches/12  inches  per  ft)  =  1300  ft/sec 

•53.35=  186.725 


X™  = 


wViprp  C’  — 

Y  ^ 

R  = 

f  1.4  ^ 

^2-Cp-g-T„ 

[y-h 

[1.4-lJ 

1300 

V2-186.725-32.2-517.S 


=  0.5210 


sin^  P,  = 


X, 


1  +  X,, 


(  Y“1  2  1 

(  1-4-1  2^ 

(  0.5210"  ^ 

1  + - 1.28" 

2 

^1+0.5210"  ^ 

1-4-1.1.28" 

t  2  J 

1  2  J 

sin"  =  0.8649 


sinP„  =0.9300  ^  P„  =  sin-‘(0.9300)  =  68.438" 
=M  •  cos  P,,  =  1 .28  •  cos(68.438")  =  0.4704 
=0.455;  R.^  =  1.0687;  =  0.956 

0, „=  1.52;  f„(a,./qj=  0.8553 


5^  =0.0045;  AR^  =1.2 

=  0.037  at  tip;  0.056  at  mean;  0.08  at  hub 
D3__^  =  0.52;  R33  =1.0259;  (^,3  =  1.0371 
3m  =1-25;C(c73,  /a,„)  =0.8 
83  =0.003;  AR3  =1.2 

(t/c)3_^^3  =0.07  at  tip;  0.06  at  mean;  0.05  at  hub 

Applying  equations  from  Appendix  A  yields  the  followings  results. 


95 


1 


Inlet  Conditions 


0  4704 
cosO 


X,  = 


1.4-1 


0.4704" 


1  4-1 

1  +  -!^^ -0.4704" 


=  0.2059 


=  (0.2059)  •  cosO  =  0.2059 
X,j,  =(0.2059)  •tan0  =  0 

=  0  +  (0.2509)  •  tan(68.438)  =  0.5210 
14.69 


P..  =; 


.  144  =  0.07658  Ibm/ft" 


53.35-517.J 

V,  =  V2- 32.2 -186.725 -517.8  =  2495.3139  ft/sec 


4>  „  =  0.2059  -  ( 1  -  0.2059" )  >  ‘^‘  =  0. 1 847 


A,  = 


17.09 


0.07658-2495.3139 


1 


0.1 847 -cosO 


-  =  0.4841  ft"  or  69.7144  in" 


0.5210-2495.3139  ^ 

r  = - =  0.4583  ft  or  5.5  in 

2836.33 

r,h  =  .^5.5" -69-71^  =  2.8389  in  or  0.2366  ft 


=  2-838 V  -0.5162 
""  /5.5 


5.5  +  2.8389 

r,„  = - : - =4.1695  in  or  0.3475  ft 


X„,  =- 


2 

5.5 

4.1695 


-0  =  0 


X^,  =X,  =0.2058 

Zlm  Zlt 


a,„  =  tan 


5.5 


V 


4.1695 
4.1695 


--  tanO 


=  0 


5.5 


-0.5210=0.3949 


Pi„  =  tan- 


0.3949 -0 
0.2058 


=  62.4685“ 


X,  =  =0.2058 

cosO 
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=  0.4704 


2.  Rotor  Conditions  at  Mean  Line 


X  =  0.4663  •  0.42 1 3  =  0. 1 968 

=0.4221-0.1968tan47.9895  =0.2036 
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a,„  =  tan 


0.2036 

0.1968 


=  45.9677° 


r.  = 


1-1 

2 


0.0  0.2036 


0.3949  0.4221 


=  0.7588 


T  =  1  +  2  •  [0.422 1-  0.2036-  0.3949  ■  O.O]  =  1 . 1 7 1 8 


0.1968 

= - =  0.2832 

cos(45.9677) 


0.1968 

= _ _  =  0.2941 


= 


cos(47.9894) 
0.2832 


Vl.1718 


=  0.2616 


0.2941 

=  ,  =0.2717 


Vl.1718 


M,„  =• 


1.4-1 


-0.2616" 


1-0.2616" 


=  0.6060 


^W2m 


1.4-1 


-0.2717" 


1-0.2717" 


=  0.6312 


Stator  Conditions  at  Mean  Line 


a  =10=0" 


^  cos45.9677 

D.  =  1-1.0371 - + 

cosO 

( tan45.9677  - 1 .0259- 1 .0371  •  tanO)  •  cos45.9677 
(1  +  1.0259) -1.25 


=  0.5631 


X„  =1.0371  0.1968  =0.2041 

Z3m 


X„,  =  0.2041  •tan0=0 

63m 


X,„  =0.42211.0259  =  0.4330 

U3m 


0.4330-0 


X,  = 


V 

0.2041 


0.2041 
=  0.2041 


=  64.7615° 


X,,.,  = 


cosO 

0.2041 
cos64.7615 


■  =  0.4787 
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=  0.2616  •  (1-0.2616" =  0.2191 


A  /  0.1847-cosO  Vl.1718 

/A,  0.2191  00845.9677  1.7218 


=  0.7627-69.7144=53.1725  in"  or  0.3693  ft" 
r  =  4.1694- 1.0687  =  4.4559  in  or  0.3713  ft 

2ni 

53.1725 


H  = 


2-71-4.4559 

^  1.8992 

y  _  2-4.4559  _  Q 

1.8992 


=  1.8992  in  or  0.1583  ft 


1  +  - 


2-4.4559 


-4.4559  =  5.4055  in  or  0.4505  ft 


1  +  0.6487 

\  y 

r  ,  =  0.6487  -  5.4055  =  3.5063  in  or  0.2922  ft 


Stator  Performance 


Initially  set  cOj  and  cbsp^^  =  0 

cos"  45.9677 

(0  =21.25 - ; - 

'  cos  0 

cb^  =0  +  0  +  0.02546 
P, 


[0.005+0.16-0.5631'*]  = 


=  1.7218 -0.02546- (1.7218-1.3435)  =  1.7121 


=  x  =  1.1718 


=  1.1718-0.2041"  =  1.1302 


/P,  = 


1.1302 

1.1718 


=  1.5084 


=  0.1886  (1-0.1886")“  =  0.1722 


A,  /  _  0.1847-cos0  Vl.1718 
/K  0.1722- cosO  1.7121 

A,  =0.6782-69.7144  =47.2778  in^  or  0.3233  ff 

n  =1.0259  -4.4559  =4.5713  in  or  0.3809  ft 


47  '?11R 

H  =  =1.6460  in  or  0.1372  ft 


2-71-4.5713 


0.02546 
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1- 


1.6460 


^ - 2  ■4.57 13  =0.6949 


!  +  ■ 


1.6460 

2-4.5713 


1  +  0.6949 


•4.5713  =  5.3943  in  or  0.4495  ft 


n,  =  0.6949-  5.3943  =  3.7483  in  or  0.3124  ft 


Rotor  Hub  Calculations 


4  4559 

=  0.2036  ■  =0.2587 

3.5063 


=  0.4221  =0.3321 

4.4559 


a,.  =  tan 


0.2587 

0.1968 


=  52.7381” 


p2h  =  tan^ 


0.3321-0.2587 

0.1968 


=  20.4569" 


X,=  =0.3251 

COS52.7381 


Xw2h  =  =0.2101 

COS20.4569 


0.3251 

Y,,  =  .  =0.3003 


•v/l.1718 


0.2101 

=  .  =0.1940 


-s/l.1718 


= 


1.4-1 


•0.3003^ 


1-0.3003" 


=  0.7039 


M 


1.4-1 


•0.1940" 


1-0.1940" 


=  0.4423 


Rotor  Tip  Calculations 


4  4559 

X„,  =0.2036  =0.1678 

5.4055 


X,„  =0.4221-^^^=0.5120 
4.4559 
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a,,  =  tan 


0.1678 

0.1968 


=  40.4526° 


p2,  =  tan  ' 


0.5120-0.1678 


V 


0.1968 


=  60.2402° 


X,.  =  =  0.2586 


Xw.= 


COS40.4526 
0.1968 


cos60.2402 


=  0.3965 


0.2586 

V ,  =  ,  =  0.2389 


Vl.1718 


Y„.,=^££fL=  0.3663 


Vl.1718 


M,  ,=■ 


1.4-1 


•0.2389" 


1-0.2389" 


•  =  0.5502 


Mw2t  ^ 


1.4-1 


•0.3663" 


1-0.3663" 


•  =  0.8802 


Stator  Hub  Calculations 


45713 

=0.0-  =0.0 
3.7483 


o  7^00 

=0.4330  •  =0.3550 

4.5713 


a,.  =tan 


0.0 

0.2041 


=  0.0° 


p3h  =  tan^ 


0.3550-0.0 
^  0.2041  ^ 


=  60.1058° 


X,,  = 


0.2041 


X„„.,  = 


cosO 
0.2041 


=  0.2041 


Y,,  = 


cos60.1058 
0.2041 


=  0.4095 


Y„,.  = 


Vl.1718 

0.4095 


=  0.1886 


Vl.1718 


=  0.3783 
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=  0.4293 


M,,  =' 


•0.1886" 


1-0.1886" 


Mw3h 


•0.3783" 


1-0.3783" 


=  0.9139 


Stator  Tip  Calculations 


4  5713 

=0.0-  =0.0 
5.3943 


5  3943 

X,„  =0.4330-  =0.5110 

4.5713 


a„  =  tan  ‘  — ^ -  =0.0° 

0.2041 


„  ,(0.5110-0.0  1 

P,,  =  tan  -  =  68.2250° 

0.2041 


X3,  =0-2041  ^0  2041 
cosO 


Xw3,=  =0.5502 

COS68.2250 


0.2041 

Y3,  =  ,  =0.1886 

■Jl.1718 


0.5502 

Y„3,  =  .  =  0.5083 

•Jl.1718 


— 0.1886" 

M3,  =  1  - =  0.4293 

1-0.1886" 


^W3t 


•0.5083" 


1-0.5083" 


=  1.3197 


Rotor  Blade  Geometry 


H,  =5.4055-3.5063=  1.8992  in 


1.8992 

C,  = - =  1.5827  in 

1.2 


2  •7t-4  559 

Z  =  -1.52=26.8884  blades 

1.5827 


Choose  Z2  to  match  transonic  compressor  rotor.  Therefore,  Z2  =  22. 
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1.8992 


AR  = 

2 -Tt- 4.559 -1.52 


=  0.9818 


22 


1  8992 

,  =  =1.9343  in 

0.9818 


1  9343 

=  =  l  2726in 

1.52 

Stator  Blade  Geometry 

H,  =  5.3943-3.7483  =  1 .6460  in 
,  1.6460 

C,  = - =  1.3717  in 

1.2 

2-7t-4.5713 


z;  = 


•1.25=  26.1740  blades 


1.3717 

Choose  Z3  to  match  transonic  compressor  rotor.  Therefore,  Z3  =  27. 
1.6460 


AR  =  - 

2 -Tt- 4.5713 -1.25 


=  1.2379 


27 


Rev3 


1.2379 
1  3297 

^  =1.0638  in 

1.25 

The  design  calculations  have  to  iterated  until  the  losses  converge  to  a  chosen  criteria.  A  criterion 
of  0.0001  was  chosen  for  both  the  hand  and  computer  calculations. 

5.  Iteration  #2 


Only  the  results  will  be  shown  for  all  subsequent  iterations. 


ROTOR 

STATOR 

56.560435987322300 

27.347157421794000 

Cl 

0.586010932091710 

Cl 

1.470018487969480 

Coi 

0.014452454402386 

Coi 

0.071059015271782 

C0SFTC2 

0.028049370920909 

WSFTC3 

0.061233413235418 

y* 

0.908770674133602 

y* 

1.978417070245130 

Pte/Pti 

0.999746596870648 

Pte/Pti 

1 .005788072708940 

C0S2 

0.000470350655165 

WS3 

-0.026347850522881 

C0p2 

0.025691775310535 

C0P3 

0.025461335476916 

cor2 

0.054211496886609 

cor3 

0.086694748712334 

Ti/T,i 

0.957618541339988 

Pt3/Pt1 

1 .666977507638590 

Pi/Pti 

0.859357579736388 

Tt3/T,i 

1.171843486972510 

Tri/T,i 

1.15597584087194 

Ts/Tti 

1.130182173979650 
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PR1/Pt1 

1.660812680211420 

Tei/T,i 

1.178143079024150 

PE1/Pt1 

1.774978623835910 

PE2/Pt1 

1.731530543151750 

Pt2/Pt1 

1.699341511605150 

T2/T,i 

1.091667880043820 

P2/Pt1 

1.326031657810530 

^2m 

0.219098107941022 

A2/A1 

0.772787543503894 

A2  (in) 

53.874387229785800 

r2m  (in) 

4.455879688743610 

h2  (in) 

1.924283336313540 

fht2 

0.644836516957841 

ra  (in) 

5.418021356900380 

r2h  (in) 

3.493738020586840 

1.468607319851730 


0.172237922031280 


0.696539265527826 


48.558787505362500 


4.571286972682070 


1 .690633536728720 


0.687879413456575 


5.416603741046430 


3.725970204317710 


Table  Cl.  Stage  Performance  Results  (Iteration  #2) 


HUB 

TIP 

)^2h 

0.259634095360950 

Y02t 

0.167421545736055 

Xu2h 

0.330933976012678 

Xu2t 

0.513206010065825 

0C2h 

52.836925539652200 

0C2t 

40.387101220022500 

P2h 

19.914379724264200 

Pa 

60.352859069903200 

X2h 

0.325797182731520 

Xa 

0.258387141279412 

)^2h 

0.209326572089755 

XiW2t 

0.397870376597357 

Y2h 

0.300962592479341 

Ya 

0.238691026272198 

Yw2h 

0.193370204379153 

Ywa 

0.367541852288360 

M2h 

0.705691535783612 

0.549615692085776 

Mw2h 

0.440706882771685 

0.883701501259568 

Mz2h 

0.426298128621960 

0.418633587082212 

Table  C2.  Rotor  Results  (Iteration  #2) 


HUB 

TIP 

^3h 

0.000000000000000 

Xest 

0.000000000000000 

Xu3h 

0.352931481111029 

Xu3t 

0.513071730606907 

ttsh 

0.000000000000000 

0.000000000000000 

CO 

CO. 

59.957896644744600 

P3t 

68.306267419874700 

Xsh 

0.204111031041579 

Xa 

0.204111031041579 

^3h 

0.407703254036660 

^3t 

0.552181051595238 

Ysh 

0.188552229153336 

Ya 

0.188552229153336 

Yw3h 

0.376625197518214 

Ywa 

0.510089864537044 
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Msh 

0.429316175796641 

IQ 

3t 

0.429316175796641 

Mw3h 

0.909100497709792 

li 

W3t 

1.326087214551410 

bh 

0.429316175796641 

li 

Z3t 

0.429316175796641 

Tabfe  C3.  Stator  Results  (Iteration  #2) 


ROTOR 

STATOR 

H2 

1.924283336313540 

Hs 

1 .690633536728720 

C2 

1.603569446927950 

Cs' 

1.408861280607270 

Zg' 

26.538058032755800 

Z3' 

25.483562094627500 

Zg 

22.000000000000000 

Z3 

27.000000000000000 

AR' 

0.994797734160298 

AR' 

1.271407814954980 

C2 

1 .934346320096740 

C3 

1 .329733478780440 

S2 

1.272596263221540 

S3 

1 .063786783024350 

Table  C4.  Blade  Geometry  Results  (Iteration  #2) 


Compare  (b^  (iteration  #2)  and  (b^  (iteration  #1). 


J  C0t2^  ^ 

1  =  0.05421  -0.02569  =  0.02852 

1 

—  0)2^3 

1  =  0.08669-0.02546  =  0.06123 

6.  Iteration  #3 


ROTOR 

STATOR 

56.560435987322300 

27.347157421794000 

Cl 

0.586010932091710 

Cl 

1.470018487969480 

Coi 

0.014264158667470 

Coi 

0.069184495778439 

raSFTC2 

0.027683925940807 

WSFTC3 

0.059618090727575 

y* 

0.908770674133602 

y* 

1.978417070245130 

Pte/Pti 

0.999746596870648 

Pte/Pti 

1 .005788072708940 

<bS2 

0.000470350655165 

OS3 

-0.026347850522881 

C0p2 

0.025691775310535 

«>P3 

0.025461335476916 

cor2 

0.053846051906507 

COt3 

0.085079426204491 

Ti/T,i 

0.957618541339988 

Pt3/Pt1 

1.667862594069170 

Pi/Pti 

0.859357579736388 

Tt3/T,i 

1.171843486972510 

Tri/T,i 

1.15597584087194 

Ts/Tti 

1.130182173979650 

PR1/Pt1 

1.660812680211420 

P3/Pt1 

1.469387081069080 

Tei/T,i 

1.178143079024150 

^3m 

0.172237922031280 

PE1/Pt1 

1.774978623835910 

A3/A1 

0.696169632289165 

PE2/Pt1 

1.731823430894990 

A3  (in) 

48.533018761 546700 

Pt2/Pt1 

1 .699628954585760 

r3m  (in) 

4.571286972682070 

T2/T,i 

1.091667880043820 

h3  (in) 

1.689736366417590 

P2/Pt1 

1.326255955569100 

fht3 

0.688019209608545 

^2m 

0.219098107941022 

r3t  (in) 

5.416155155890860 

A2/A1 

0.772656848887118 

r3h  (in) 

3.726418789473280 
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As  (in) 

53.865275938523100 

r2m  (in) 

4.455879688743610 

hs  (in) 

1.923957899554990 

fht2 

0.644885917495964 

ra  (in) 

5.417858638521100 

r2h  (in) 

3.493900738966110 

Table  C5.  Stage  Performance  Results  (Iteration  #3) 


HUB 

TIP 

^2h 

0.259622003649605 

Xe2t 

0.167426574025716 

Xu2h 

0.330949389028739 

Xu2t 

0.513190597049764 

a2h 

52.835640920353900 

aa 

40.387950481631600 

P2h 

19.921457676992400 

Pa 

60.351402895843300 

Xsh 

0.325787546694057 

Xa 

0.258390399380250 

X«2h 

0.209335942224263 

XA/2t 

0.397852611427487 

Ysh 

0.300953690969537 

Ya 

0.238694036017300 

Yw2h 

0.193378860255026 

Ywa 

0.367525441306733 

0.705668584925989 

0.549623041103662 

Mw2h 

0.440727376594600 

0.883655881142935 

Mz2h 

0.426296872834759 

0.418633906000314 

Table  C6.  Rotor  Results  (Iteration  #3) 


HUB 

TIP 

^3h 

0.000000000000000 

Xest 

0.000000000000000 

Xu3h 

0.352973972009957 

Xu3t 

0.513029239707979 

OCsh 

0.000000000000000 

«3t 

0.000000000000000 

P3h 

59.960885862080200 

P3t 

68.304637549438200 

Xsh 

0.204111031041579 

Xa 

0.204111031041579 

X«3h 

0.407740037167485 

Xwa 

0.552141570422119 

Ysh 

0.188552229153336 

Ya 

0.188552229153336 

Yw3h 

0.376659176775861 

Ywa 

0.510053392901174 

Mah 

0.429316175796641 

IH 

0.429316175796641 

0.909196076420228 

1 .325959056752780 

Mz3h 

0.429316175796641 

Mz3t 

0.429316175796641 

Table  C7.  Stator  Results  (Iteration  #3) 


ROTOR 

STATOR 

Hs 

1 .923957899554990 

Ha 

1.689736366417590 

Cs' 

1.603298249629160 

Ca' 

1.408113638681320 
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Zg' 

26.542546935338600 

Zb' 

25.497092664122100 

Zg 

22.000000000000000 

Zb 

27.000000000000000 

AR' 

0.994629492953864 

AR' 

1.270733115606990 

Cg 

1 .934346320096740 

Cs 

1 .329733478780440 

Sg 

1.272596263221540 

S3 

1 .063786783024350 

Table  C8.  Blade  Geometry  Results  (Iteration  #3) 


Compare  (b^  (iteration  #3)  and  d)^  (iteration  #2). 
-®T2„„,„.|=  0.05385-0.05421  =  0.0003654 

=  0.08508-0.08669  =  0.001615 

7.  Iteration  #4 


ROTOR 

STATOR 

56.560435987322300 

27.347157421794000 

Cl 

0.586010932091710 

Cl 

1.470018487969480 

Coi 

0.014266571444569 

Coi 

0.069221229482487 

C0SFTC2 

0.027688608666518 

WSETC3 

0.059649745121757 

y* 

0.908770674133602 

y* 

1.978417070245130 

Pte/Pti 

0.999746596870648 

Pte/Pti 

1 .005788072708940 

ras2 

0.000470350655165 

ras3 

-0.026347850522881 

C0p2 

0.025691775310535 

K>P3 

0.025461335476916 

a>r2 

0.053850734632218 

cor3 

0.085111080598674 

Ti/Tn 

0.957618541339988 

Pt3/Pt1 

1.667847160812340 

Pi/Pti 

0.859357579736388 

Tt3/T,i 

1.171843486972510 

Tri/T,i 

1.15597584087194 

T3/Tti 

1.130182173979650 

PR1/Pt1 

1.660812680211420 

P3/Pt1 

1 .469373484368560 

Tei/T,i 

1.178143079024150 

^3m 

0.172237922031280 

PE1/Pt1 

1.774978623835910 

A3/A1 

0.696176074225205 

PE2/Pt1 

1.731819677900590 

A3  (in) 

48.533467856980100 

Pt2/Pt1 

1.699625271359260 

r3m  (in) 

4.571286972682070 

Tg/T,i 

1.091667880043820 

h3  (in) 

1 .689752002224000 

P2/Pt1 

1.326253081470570 

fht3 

0.688016773055056 

^2m 

0.219098107941022 

rst  (in) 

5.416162973794070 

A2/A1 

0.772658523297490 

rsh  (in) 

3.726410971570070 

Ag  (in) 

53.865392668966700 

rgm  (in) 

4.455879688743610 

h2  (in) 

1 .923962068928580 

fht2 

0.644885284576048 

ta  (in) 

5.417860723207900 

tgh  (in) 

3.493898654279320 

Table  C9.  Stage  Performance  Results  (Iteration  #4) 
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HUB 

TIP 

^2h 

0.259622158556950 

Xe2t 

0.167426509603247 

Xu2h 

0.330949191562962 

Xu2t 

0.513190794515540 

a2h 

52.835657378138000 

tta 

40.387939601023500 

P2h 

19.921367002758800 

Pa 

60.351421552757200 

X2h 

0.325787670140640 

X2t 

0.258390357637089 

X«2h 

0.209335822159843 

XA/2t 

0.397852839028210 

Y2h 

0.300953805006140 

Y2t 

0.238693997456101 

Yw2h 

0.193378749342774 

Ywa 

0.367525651558102 

M2h 

0.705668878946494 

0.549622946947096 

Mw2h 

0.440727113995883 

0.883656465604799 

Mz2h 

0.426296888922253 

0.418633901914269 

Table  CIO. 

Rotor  Results  (Iteration  #4) 

HUB 

TIP 

^3h 

0.000000000000000 

Xest 

0.000000000000000 

Xu3h 

0.352973231482253 

Xu3t 

0.513029980235682 

ttSh 

0.000000000000000 

0.000000000000000 

CO 

CO. 

59.960833770877400 

Pat 

68.304665956674600 

Xsh 

0.204111031041579 

Xst 

0.204111031041579 

X«3h 

0.407739396104768 

X/V3t 

0.552142258492756 

Ysh 

0.188552229153336 

Yat 

0.188552229153336 

Yw3h 

0.376658584579527 

Yw3t 

0.510054028522148 

0.429316175796641 

0.429316175796641 

Mw3h 

0.909194410621443 

1.325961290180600 

Mz3h 

0.429316175796641 

0.429316175796641 

Tabfe  Cl  1.  Stator  Results  (Iteration  #4) 


ROTOR 

STATOR 

Hs 

1 .923962068928580 

Ha 

1 .689752002224000 

C2' 

1.603301724107150 

Ca' 

1.408126668520000 

Z2' 

26.542489415600500 

Za' 

25.496856731509300 

Za 

22.000000000000000 

Za 

27.000000000000000 

AR' 

0.994631648397049 

AR' 

1 .270744874208740 

C2 

1 .934346320096740 

Ca 

1 .329733478780440 

S2 

1.272596263221540 

Sa 

1 .063786783024350 

Table  C12.  Blade  Geometry  Results  (Iteration  #4) 


Compare  cb^  (iteration  #4)  and  d)^  (iteration  #3). 

=  0.05385 - 0.05385  =  0.000004683 
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=  0.085 1 1  - 0.08508  =  0.00003 165 

Both  the  rotor  and  stator  losses  satisfy  the  criterion  after  the  fourth  iteration, 
n^,  =  —  =  1.6678 

P,. 


Tc=tp=  1.1718 


lie  =- 


1.41 

1.6678~-1 
1.1718-1 


=  0.9158  or  91.58% 


C2.  CODE  RESULTS 


HUB 

MEAN 

TIP 

Plh 

52.195176449181800 

Plm 

62.343578122316700 

Pit 

68.410458367312200 

aih 

0.000000000000000 

0tlm 

0.000000000000000 

ait 

0.000000000000000 

Xih 

0.205867575543143 

Xim 

0.205867575543143 

Xit 

0.205867575543143 

Xuih 

0.265557741712975 

Xuim 

0.393264467850964 

Xuit 

0.520971193988953 

rih 

2.794827442941060 

Um 

4.138860046006270 

rit 

5.482892649071480 

^1h 

0.000000000000000 

Xeim 

0.000000000000000 

Xen 

0.000000000000000 

0.205867575543143 

0.205867575543143 

)tit 

0.205867575543143 

Mzih 

0.470410144688393 

Mzim 

0.470410144688393 

0.470410144688393 

^1h 

0.336009483264546 

^1m 

0.443890076859253 

Xwit 

0.560171798313952 

0.470410144688393 

Hll 

0.470410144688393 

0.470410144688393 

0.797720638371130 

nBi 

1.107676776722660 

1 .280000000000000 

Ai 

69.931875719909200 

fhtl 

0.509735940829400 

CO 

2836.334567415990000 

mdot 

17.090000000000000 

14.690000000000000 

T,i 

517.800000000000000 

R 

53.350000000000000 

Y 

1 .400000000000000 

Ta 

ble  C13.  Inlet  Results 

HUB 

MEAN 

TIP 

P2h 

19.153113641505600 

P2m 

47.808100398602100 

Pa 

60.322721543344400 

0C2h 

52.844869054566300 

45.871861712689900 

0C2t 

40.226189216790300 

X2h 

0.326016514230627 

X2m 

0.282758297008589 

X2t 

0.257833958315095 

Xu2h 

0.328294695162428 

Xu2m 

0.420281736792325 

Xu2t 

0.512268778422222 

r2h 

3.455094238614300 

f2m 

4.423199731166900 

f2t 

5.391305223719500 

)^2h 

0.259909266378080 

Xe2m 

0.203022938914057 

X02t 

0.166566531027495 

no 


^2h 

0.196809402219245 

Mz2h 

0.426564947028765 

)^2h 

0.208351884253337 

M2h 

0.706608604848915 

Mw2h 

0.438807084275359 

Y2h 

0.301318214114674 

Yw2h 

0.192567599892291 

D2h 

0.539065878334816 

fst2h 

0.604152503515955 

0.196809402219245)^2, 


0.421383869310173 


0.293147277076569)^21 


0.605406977153977 


0.629379184024538 


0.261337145088495  Ya 


0.270939078683619  Ywa 


0.455000000000000  Da 


0.758468045193244  rsct 


1.170653666750960 


0.640864149819104 


1.520000000000000 


0.855300000000000 


1.068700000000000 


0.956000000000000 


Table  C 1 4 .  Rotor  Results 


0.196809402219245 


0.421383869310173 


0.397798924905110 


0.548663539641709 


0.629379184024538 


0.238301019937484 


0.367662545905120 


0.390436898079840 


0.837422718264700 


HUB 

MEAN 

TIP 

Psh 

59.760648493304300 

Psm 

64.641480976070200 

Pat 

68.231904833353900 

ttSh 

0.000000000000000 

Wsm 

0.000000000000000 

ttat 

0.000000000000000 

)^h 

0.204111031041579 

)^m 

0.204111031041579 

)^3t 

0.204111031041579 

)^3h 

0.350482583299427 

)^3m 

0.431167033775246 

)^3t 

0.511851484251066 

fSh 

3.688607742179240 

l"3m 

4.537760604204120 

fat 

5.386913466229010 

)^3h 

0.000000000000000 

)^3m 

0.000000000000000 

)^3t 

0.000000000000000 

)t3h 

0.204111031041579 

)t3m 

0.204111031041579 

)^3t 

0.204111031041579 

Mz3h 

0.429542339664783 

Mz3m 

0.429542339664783 

Mz3t 

0.429542339664783 

)^3h 

0.405585199667217 

)^3m 

0.477039122093147 

)^3t 

0.551047416220125 

Mah 

0.429542339664783 

0.429542339664783 

0.429542339664783 

0.904138170145811 


0.188648024471074 


0.374858949472163 


0.636269323766171 


0.604152503515955 


1.098406365761990 


0.188648024471074  Ys, 


0.440899678567374  Ywst 


0.561 67486245301 5  Ds, 


0.7584680451 93244  rs,3t 


1.000000000000000 


0.684734916442117 


1.250000000000000 


0.800000000000000 


1.025900000000000 


1.037100000000000 


1.098406365761990 


0.188648024471074 


0.509301265734333 


0.527244214226719 


0.837422718264700 
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ROTOR 

raSFTC2 

0.015092502816179 

C0S2 

0.000415286496133 

C0p2 

0.025621919830276 

a>r2 

0.041129709142588 

Ti/T,i 

0.957618541339988 

Pi/Pti 

0.859357579736388 

TRi/Tn 

1.15465694167410 

PR1/Pt1 

1.654190026081820 

Tei/T,i 

1.176636738281170 

PE1/Pt1 

1.767048277556160 

PE2/Pt1 

1 .734357050220880 

Pt2/Pt1 

1.703686175109830 

T2/T,i 

1.090701412223770 

P2/Pt1 

1.330027617869200 

^2m 

0.218975111752094 

A2/A1 

0.769783662636916 

A2  (in) 

53.832415426741300 

STATOR 


0.052221964041667 


0.000000000000000 


0.025342625084949 


0.077564589126615 


1.674703502643810 


1.170653666750960 


1.128992353758110 


1 .475220458598720 


0.172309288336835 


0.692686802248097 


48.440887367635200 


1.674703502643810 


1.170653666750960 


1.170653666750960 


Table  Cl 6.  Stage  Performance  Results 


ROTOR 

STATOR 

H2 

1.936210985105210 

H3 

1 .698305724049770 

C2' 

1.613509154254340 

C3' 

1.415254770041480 

Z2' 

26.191679106060200 

25.192519448425300 

Z2 

22.000000000000000 

27.000000000000000 

AR' 

1.007953705186150 

AR' 

1.277177533055250 

C2 

1 .920932454678200 

C3 

1.320512345138770 

S2 

1.263771351761970 

S3 

1.056409876111010 

i*2h 

3.027273612963370 

i*3h 

-14.428782843997100 

I* 

1  2m 

6.757152202149710 

I* 

1  3m 

-1.724296488710650 

I* 

1  2t 

10.291562745234000 

I* 

1  3t 

3.241975399750170 

CM 

-42.100896728952500 

<t>*3h 

76.083582934974900 

(t>*2m 

-24.523123063313700 

3m 

26.611448616369400 

<t>*a 

-19.983188007865800 

<t'*3t 

8.263624315276530 

5*2h 

-5.757137628140870 

21 .988386237067600 

5*2m 

-2.833187565346610 

8.688685825842080 

5*a 

-1.152333894067850 

5*3t 

4.391711979088870 

Table  C17.  Blade  Geometry  Results 
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Ratio/Efficiency 


APPENDIX  D.  COMPARISON  OF  RESULTS 


Figure  Dl.  Annulus  Comparison 


n 

X 

T) 

Stage  Performance 

L 

DHand  Calc.  Results 

o  Code  Results  ^  Transonic  Compressor  I 

Figure  D2.  Stage  Performanee  Comparison 
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Inches 


2.50 


2.00 


1.50 


Figure  D3.  Solidity  Comparison 


□  Hand  Calc.  Results_ HCode  Results □  Transonic  Compressor 


Figure  D4.  Chord  and  Blade  Height  Comparison 
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1000 


900 

800 

700 


600 


VZ2h 

VZ2m 

VZ2t 

V2h 

V2m 

Velocity 

D  Hand  Calc.  Results 

aCode  Results 

□  Transonic  Compressor 

Figure  D5.  Rotor  Velocity  Comparison 


□  Hand  Calc.  Results  b  Code  Results  □  Transonic  Compressor 

Figure  D6.  Rotor  Flow  Angle  Comparison 
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Ft/sec 


rst2h  rst2m  rst2t  D2h  D2m  D2t 

Deg.  of  Reaction/Diffusion 

I  °  Hand  Calc.  Results  °  Code  Results  °  Transonic  Compressor  | 


Figure  D7.  Rotor  Degree  of  Reaction  and  Diffusion  Comparison 


VZ3h  VZ3m 

VZ3t 

V3h 

V3m 

Velocity 

1  o  Hand  Calc.  Results 

B  Code  Results 

DTransonic  Compressor 

Figure  D8.  Stator  Velocity  Comparison 
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